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FOREWORD 


This  report  is  one  of  a  series  presenting  the  results  of  research 
performed  in  the  field  of  shrouded  propeller  aerodynamics,  under  contract 
for  the  U.  S.  Army  Transportation  Research  Command,  Fort  Eustis,  Virginia. 

The  work  performed  by  Dr.  Theodorsen  and  his  associates  at  the  Republic 
Aviation  Corporation  is  presented  in  four  volumes.  The  first  two  volumes 
were  issued  during  August  1960,  under  Contract  DA  44- 177-TC-606.  This 
volume  presents  the  theoretical  approach  to  ducted  propeller  flow,  while 
Volume  IV  presents  a  discussion,  of  vertical  take-off  vehicles. 

The  report  has  been  reviewed  by  the  U.  S.  Army  Transportation  Research 
Command  and  is  considered  to  be  technically  sound.  The  report  is  published 
for  the  exchange  of  information  and  the  stimulation  of  ideas. 
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INTRODUCTION 


As  this  investigation  is  limited  to  ducted  propeller  aerodynamics 
the  following  conclusions  are,  in  general,  restricted  to  this  field  with 
reference  to  other  types  for  purposes  of  comparison  or  to  emphasize 
differences.  The  present  treatment  is  perfectly  general  and  not  intended 
to  cover  any  particular  design  or  to  give  elementary  design  data.  The 
purpose  is,  on  the  other  hand,  to  illuminate  the  fundamental  knowledge 
required  to  arrive  at  a  preliminary  layout  which,  at  a  later  stage,  does 
not  provide  undesirable  and  unexpected  surprises. 

It  is  perfectly  clear  that  such  a  considerable  departure  from  the 
theory  of  conventional  aircraft  can  be  handled  effectively  only  by  engineers 
or  scientistswell  skilled  in  the  treatment  of  complex  airflow  problems. 

It  should  be  noted  that  the  actual  performance  can  be  predicted  with 
considerable  precision  since,  with  a  few  exceptions,  no  really  unknown 
factors  ctre  involved.  Let  us  take  as  an  example  the  propeller  or  fan 
itself.  In  the  first  place,  it  is  perfectly  well  known  that  the  propeller 
efficiency  reaches  its  maximum  only  with  uniform  circumferential  and 
a  prescribed  and  predictable  radial  load  distribution.  Since  the  propeller 
is  or  may  be  located  in  a  duct,  it  is  quite  obvious  that,  in  addition,  the 
flow  pattern  of  the  duct  itself  must  be  known. 
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This  particular  problem  has  been  treated  in  detail  in  Ref.  Al, 

Chapter  IV  and  the  universal  flow  lines  are  shown  in  Figures  I  and  II  of 
this  reference.  These  are  the  exact  flow  lines  which  correspond  to  the  ideal 
case  of  uniform  downflow  velocity  at  the  plane  of  the  propeller.  Great 
confusion  has  existed  on  this  problem  due  to  a  lack  of  basic  understanding. 

By  inspection  of  the  pattern  of  the  ideal  streamlines,  any  one  of  which  may 
represent  the  wall  of  the  duct  inlet,  it  should  be  noticed  that  any  smaller 
radius  of  curvature  of  the  duct  entrance  results  unavoidably  in  an  increase 
in  the  inlet  velocity  at  and  near  the  wall.  A  standard  type  propeller  operating 
in  this  inlet  will  thus  be  unloaded  at  the  tip  producing  a  condition  quite 
detrimental  to  the  efficiency  of  the  propeller.  Conversely,  the  inlet  must  be 
designed  as  shown  in  Ref.  Al  in  order  to  reach  maximum  efficiency  with  a 
normal  type  propeller,  that  is  a  propeller  designed  to  operate  in  a  circular  duct. 

The  second  step  in  the  procedure  is  then  to  obtain  the  load  distribution 
on  a  multiblade  propeller  in  a  circular  infinitely  long  duct.  The  theory  of  the 
optimum  loading  of  propellers  has  been  given  by  Prandtl,  and  with  more 
refinement  by  Goldstein  (Ref.  A2).  For  counter  rotation  propellers  which  are 
desirable  for  heavier  loading,  the  general  theory  is  given  by  Theodorsen  (Ref. 
A3).  However,  the  actual  case  of  a  propeller  in  a  circular  duct,  single  or 
counter  rotating,  has  not  been  specifically  treated  in  any  of  the  above  refer¬ 
ences,  being  of  no  significance  at  the  time.  The  single  rotating  case  has  been 
treated  in  the  present  investigation  in  Section  B.  Both  the  loading  functions 


and  the  mass  coefficients  are  expressed  by  appropriate  formulas. 
Numerical  results  and  graphs  remain  to  be  processed.  With  such 
tables  and  graphs  available  for  multiblade  propellers  the  ideal  propeller 
can  be  designed  with  perfect  accuracy. 

It  may  be  remarked  that  although  the  duct  -streamlines  are 
independent  of  the  magnitude  of  the  flow  velocity,  the  propeller  itself 
has  to  be  designed  for  a  prescribed  velocity  and  thrust.  As  the  take¬ 
off  condition  requires  the  highest  power  and  efficiency,  the  propeller 
twist  distribution  must  be  designed  to  fit  this  case.  If  a  variable  pitch 
propeller  is  employed,  the  loss  in  efficiency  is  not  excessive  as  the  thrust 
is  being  reduced.  However,  any  decrease  in  the  angle  of  attack  of  the 
propeller  will  decrease  the  loading  at  the  tip  relative  to  that  at  the  root 
section.  This  is  particularly  true  if  the  design  employs  a  relatively  low 
lift  coefficient. 

In  summary,  to  obtain  the  most  efficient  duct-propeller  com¬ 
bination,  the  duct,  whenever  possible,  should  be  designed  with  internal 
flow-lines  so  designed  as  to  provide  a  uniform  flow  over  the  inlet  area 
at  the  location  of  the  propeller.  The  propeller  will  then  be  designed  for 
the  ideal  load  distribution.  Experimental  checking  either  on  full  scale  or 
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on  models  is  not  at  all  required,  if  and  when  the  ideal  combination 
may  be  employed. 

If  and  when  deviations  from  the  ideal  case  of  duct  design  are 
unavoidable,  it  is  still  possible  to  estimate  very  nearly  the  resulting 
radial  non- uniformity  of  the  inflow  and  correct  the  blade  angle  distri¬ 
bution  of  the  propeller  accordingly.  Only  if  drastic  deviations  from 
the  ideal  inlet  are  employed  may  it  be  desirable  to  run  a  final  ex¬ 
perimental  check  rather  than  to  obtain  the  services  of  an  organization 
sufficiently  qualified  to  obtain  the  ideal  combination  by  direct  calculations. 

As  a  further  remark  it  is  stated  that  the  complimentary  combination 
always  exists:  for  any  duct  design,  the  complimentary  propeller  is  a 
matter  of  proper  design.  This  point  must  not  be  overlooked,  since  any 
deficiency  reflects  itself  directly  in  the  weight  of  the  power  plant  and  in 
the  fuel  consumption. 

As  we  are  not,  in  particular,  referring  to  any  specific  design,  we 
shall,  nevertheless,  note  that  there  exists  somewhat  of  a  difference  between 
the  case  of  a  propeller  in  a  wing  and  the  case  of  a  free  duct  a-la-Doak.  In 
the  latter  case  the  numerical  calculation  of  the  flow  velocity  in  the  inlet 
plane  in  the  take-off  condition  is  quite  complex.  However,  it  is  possible 
to  estimate,  at  least  tc  the  first  order  of  magnitude,  the  effect  of  any 
particular  radius  of  curvature  at  the  duct  inlet.  The  velocity  at  the  propeller 
tip  will  be  somewhat  greater  than  the  average  and  a  corresponding  correction 
must  be  applied  to  the  propeller  twist  distribution.  Experimentally,  if 


theoretical  determination  is  unavailable,  one  need  only  insert  the  duct  as 
an  extension  of  a  long  cylindrical  pipe  of  the  same  diameter  as  that  of  the 
duct  exit,  apply  a  suction  at  the  other  end  and  measure  the  radial  velocity 
distribution  at  the  plane  of  the  location  of  the  intended  propeller.  With 
the  other  data  available  the  propeller  may  thus  be  fully  specified.  It 
should  be  emphasized  that  the  duct  enclosing  the  propeller  is  not  necessarily 
beneficial.  We  shall  discuss  this  matter  briefly.  (See  also  Volume  IV. ) 

Favoring  the  duct  arrangement  are  the  following  facts: 

1.  The  propeller  diameter  is  reduced.  The  theoretical 
limit  of  such  reduction  as  shown  in  the  first  report 
(Ref..Al)  is  in  the  order  of  30%  based  on  the  diameter 
of  the  free  propeller,  but  actually,  depending  on  the 
case,  only  about  10  to  15%. 

2.  The  RPM  of  the  propeller  and  the  related  shafting  and 
gearing  is  (correspondingly)  higher  resulting  in  a  saving 
of  weight. 

The  adverse  factors  are  as  follows: 

1.  Skin  friction  of  the  duct  walls  has  been  added  to  that  of 
the  assembly.  Based  on  the  area  of  the  duct,  the 
solidity  of  the  propeller  and  the  flow  velocities,  it  will 
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be  found  that  the  duct  loss  generally  outweighs  the 
reduction  in  propeller  skin  friction  loss.  Again  there 
is  a  limit:  If  the  propeller  is  heavily  loaded  a  balance 
may  result.  No  credit  should  be  given  to  the  duct  as 
a  lifting  surface  since  the  critical  situation  exists  at 
take  off  and  in  transition,  in  which  cases  any  potential 
lifting  capacity  is  destroyed  by  stalling. 

2.  Adverse  is,  of  course,  also  the  effect  of  the  weight 

of  the  duct.  This  must  be  balanced  against  the  decrease 
in  weight  of  the  shafting  and  gearing.  Adverse  is  also 
the  external  drag  of  the  duct. 

The  effects  of  transition  or  forward  speed  will  be  discussed 
next.  It  is  quite  essential  to  be  aware  of  these  effects  in  the  early 
design  stage  to  avoid  later  difficulties.  The  most  essential  require¬ 
ment  is  an  understanding  of  the  nature  and  magnitude  of  the  various 
effects. 


Pitch-up  moments  are  the  first  in  line  to  be  considered.  Large 
pitch-up  moment  represents  an  inherent  deficiency  of  the  propeller-in¬ 
wing  arrangement.  It  has  been  shown  in  the  previous  investigation 
(Ref.  Ai,  Chapter  II)  that  for  a  propeller  inserted  in  a  surface  of  area 
A  there  exists  an  associated  pitch-up  moment  equal  to  the  full  momentum 
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drag  times  a  distance  equal  to  one  half  of  the  "mean  radius"  of  the 
area  A  and,  as  may  be  noted,  fairly  independent  of  the  shape.  It 
should  be  pointed  out  that  this  moment  cannot  be  eliminated  by  baffles 
or  any  other  means  as  has  sometimes  been  attempted. 

In  the  case  of  shrouded  propeller  a-la-Doak,  the  case  is  less 
serious.  The  adverse  pitching  moment  is  reduced  with  forward  tilt  of 
the  propeller  axis  as  has  been  shown  in  Ref.  A.  (Page  106).  Also  the 
numerical  value  of  the  moment  is  smaller.  The  momentum  of  the 
inlet  air-column  attacks  with  an  arm  of  approximately  one  half  the  radius 
of  the  bellmouth  at  zero  tilt  angle.  As  the  tilt  angle  is  increased  the  arm 
is  reduced  to  offset  the  gradual  increase  in  the  absolute  magnitude  of  the 
momentum  drag  with  forward  speed.  Also,  the  propeller  itself  has  a 
slightly  favorable  pitching  moment.  This  is  caused  by  the  fact  that  the 
air  enters  the  propeller  disk  with  more  downward  velocity  at  the  front 
edge  and  correspondingly  less  velocity  at  the  rear.  This  contribution 
to  cancel  part  of  the  pitch-up  moment  is  relatively  more  effective  in  the 
case  of  the  separate  duct. 

Finally,  the  pitch-up  moment  of  a  free  propeller  may  be  mentioned. 
From  the  theory  of  a  lifting  surface  it  can  be  concluded  that  the  lift  force 
is  as  always  concentrated  near  the  front  quarter  chord.  This  fact  has  been 
confirmed  by  tests  (See  Ref. A,  Page  106  and  Figures  1 6 1  and  165).  The 
magnitude  of  the  pitch-up  moment  is  quite  considerable  but  with  the  lift 
force  maintained  constant  the  value  of  the  moment  does  not  change 
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appreciably  with  forward  speed  so  in  this  hypothetical  case  the  effect 
may  be  eliminated  by  choosing  the  proper  location  of  the  hypothetical 
free  lift  propeller  with  respect  to  the  center  of  gravity  of  the  aircraft. 

To  summarize:  The  propeller-in-wing  type  suffers  from  an 
inherent  adverse  pitching  moment  of  considerable  magnitude.  The 
numerical  values  may  be  obtained  with  considerable  accuracy  by  use 
of  the  methods  given  in  the  previous  report  (Ref.  Al).  The  propeller- 
in-duct  case  is  more  favorable  as  the  acting  moment  arm  is  smaller 
and  is  gradually  decreased  in  transition  to  forward  speed.  The  hypo¬ 
thetical  case  of  a  free  propeller  is  essentially  free  of  the  defect. 

Transition  effects  on  propeller  or  fan  efficiency  is  a  serious 
matter.  The  problem  is  closely  related  to  the  problem  of  propeller 
operating  life  as  function  of  vibratory  stresses.  A  free  propeller  is 
thus  out  of  the  question,  and  a  helicopter  design  is  necessary.  The 
basic  intent  of  the  duct  design  is,  of  course,  to  force  the  flow  to 
become  more  realigned  with  the  propeller  axis. 

As  contrasted  with  the  case  of  a  free  propeller,  the  opposite 
case  is  a  propeller  with  a  long  entrance  duct.  The  propeller  would 
then,  if  designed  according  tp  theory,  reach  full  efficiency  and  would 
not  be  exposed  to  vibratory  stresses  caused  by  the  unsymmetric  flow 
pattern.  As  the  duct  may  be  shortened,  the  efficiency  of  the  lifting  fan 
will  gradually  decrease  as  a  function  of  the  nonsymmetry  of  the  flow 


pattern.  The  center  of  lift  of  the  fan  will  gradually  move  from  the 
axis  toward  the  rear  and  laterally  into  the  advancing  quadrant.  Anyone 
skilled  in  the  art  can  calculate  with  adequate  precision  the  effect  of  a 
given  travel  on  the  efficiency  of  the  propeller,  Also,  the  corresponding 
one-per-turn  vibratory  stresses  can  be  estimated  with  adequate  accuracy. 

In  regard  to  the  required  length  of  a  duct  it  is  again  clear  that  the 
case  of  the  fan  in  the  wing  is  again  in  the  most  unfavorable  position.  The 
fan-in-fuselage  and  the  separate  ducts  are  more  favorable.  Another 
parameter  of  equal  importance  is  the  ratio  of  forward  speed  to  that  of  the 
fan  tip  velocity.  A  low  value  of  the  aircraft  forward  speed  to  the  fan  inlet 
velocity  is,  of  course,  beneficial.  An  example  of  such  favorable  combi¬ 
nation  is  thus  the  GE  high  pressure  fan  type  when  installed  in  a  fuselage. 

To  prevent  excessive  vibratory  stresses  and  to  improve  the  fan 
efficiency  it  is  desirable  to  employ  inlet  vanes.  These  cannot  be  treated 
in  a  general  case  and  present  quite  a  difficult  construction  problem  if  they 
are  to  give  high  efficiency  in  the  entire  transition  range.  An  example  of  a 
simple  deflection  grid  is  treated  in  Sections  G,  H,  and  I  of  the  present 
report. 

Another  effect  apparent  in  the  transition  to  forward  flight  is  an 
interference  effect  between  the  fan-in-wing  and  the  wing  itself.  In  Ref.Al, 
Chapter  III  is  presented  the  classical  solution  of  a  wing  with  a  sink  on  the 
upper  surface  and  a  line  jet  issuing  from  the  lower  side.  These  results 


are  exact  and  show  that  the  circulation  is  increased  by  the  action  of 
the  propeller  in  an  amount  given  by  the  expression  III- 18  on  Page  25 
of  reference  report.  The  direct  contribution  of  the  propeller  is  given 
by  the  standard  expression  III- 19  on  Page  26.  Finally,  the  classical 
expression  for  the  moment  arm  of  the  pitch-up  moment  in  terms  of  the 
half  chord  of  the  wing  is  given  for  the  two-dimensional  case  by  the 
expression  III-21  on  Page  27. 

We  shall  next  consider  another  very  important  interference 
effect.  In  the  theoretical  treatment,  Ref.Al,  the  jet  was  treated  as  a 
line  jet  issuing  from  the  lower  surface  in  order  to  obtain  the  essential 
facts  of  a  complex  problem.  However,  the  jet  is  of  considerable 
dimensions  and  in  the  present  report  the  displacement  jet  on  the  wing 
proper  has  been  investigated.  The  two-dimensional  pattern  of  the  flow 
around  a  cylinder  and  the  pressure  distribution  on  a  flat  plate  perpen¬ 
dicular  to  the  cylinder  is  also  known.  By  reference  to  Section  E  of  this 
work  it  may  be  seen  that  the  loss  of  lift  on  the  area  adjacent  to  the  jet 

approaches  the  value  C  «  -  1  based  on  the  area  of  the  jet  itself.  This 

.Li 

is  a  theoretical  fact:  There  is  a  lift  reduction  caused  by  the  displacement 
of  the  jet  which  has  to  be  or  may  be  recovered  by  a  noticeable  change  in 
the  angle  of  attack  of  the  wing.  If  the  wing  area  is  ten  times  the  jet  cross 
section,  the  angle  of  attack  of  the  main  wing  must  be  increased  by  about 
one  degree.  This  gives  an  indication  of  the  magnitude  of  this  loss  of  lift. 
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A  more  serious  effect  is,  however,  that  of  the  negative  pressure 
region  existing  behind  the  jet  as  a  result  of  the  breakdown  of  the  air  flow 
around  the  jet.  An  investigation  has  been  conducted  under  Section  F  to 
clarify  the  situation,  based  on  the  theoretical  study  in  Section  E  and  on  a 
series  of  experiments  conducted  by  NASA  on  cylinders  and  jets  at  right 
angles  to  a  flat  plate  and  with  the  airstream  perpendicular  to  the  cylinder 
or  jet.  Incidentally,  the  data  confirm  closely  to  the  theoretical  data  for 
the  front  area  ahead  of  the  cylinder.  Behind  the  jet  exists,  however,  a 
large  negative  pressure  region  resulting  in  a  negative  lift  coefficient  of 
CT  =  -3  and  an  adverse  moment  coefficient  of  about  C,,  :  6  based  on  the 
area  and  the  radius  of  the  jet,  respectively.  These  values  are,  however, 
obtained  for  low  Reynolds  number  below  the  critical  value  of  the  drag 
coefficient  of  a  cylinder. 

There  is  little  doubt  that  the  negative  area  behind  the  jet  is,  to  a 
large  extent,  caused  by  the  pumping  action  of  the  jet  which  blows  the 
passing  main  airstream  downward  due  to  the  mixing  at  the  intersurface. 
The  normal  negative  region  observed  behind  a  solid  cylinder  is  thus 
intensified  by  the  downward  momentum  imposed  on  the  airstream  by  the 
jet.  The  presence  of  any  extended  surface  behind  the  jet  will,  therefore, 
contribute  to  the  magnitude  of  pitch-up  moment.  On  the  whole  it  must 
therefore  be  stated  that  the  fan-in-wing  design  is  not  justified  from 
aerodynamic  reasons  only. 
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RECOMMENDATION  ON  FURTHER  RESEARCH 


There  appears  to  be  at  least  two  distinct  problems  which  might 
advantageously  be  subject  to  further  research.. 

The  first  field  is  the  calculation  of  propeller  load  distribution 
for  a  multiblade  propeller  in  a  circular  duct,  based  on  the  formulas  of 
Section  B  of  the  present  volume. 

The  second  problem  is  the  basic  study  over  a  wide  range  of 
Reynolds  number  of  the  obstructive  effect  and  the  "pumping  action"  of 
a  jet  on  the  flow  pattern  and  pressure  distribution  on  the  lower  surface 
of  an  airfoil  or  more  generally  on  any  extended  plane  surface  parallel 
to  the  airstream.  This  work  would  be  conducted  and  would  represent 
an  extension  to  the  study  reported  in  Section  F  and  would  serve  as  a 
master  reference  for  special  cases  and  might  possibly  be  of  such 
general  value  as  to  obviate  investigations  of  individual  cases  as 


proposed  above. 


SECTION  A 


GENERAL  METHODS  FOR  THE  DESIGN  OF  VERTICAL  TAKE-OFF 

AIRCRAFT 

In  the  following  is  given  a  condensed  summary  of  the  design 
information  that  has  been  produced  by  the  work  of  the  present  and  the 
immediately  preceeding  investigation  reported  in  Ref.  Al.  The  work 
covers  the  aerodynamic  aspects  of  the  problem  and  is  presented  under 
separate  headings  indicated  below.  The  purpose  of  the  presentation  is 
to  provide  a  basis  for  rational  design  methods  based  on  a  theoretical 
analysis  of  each  problem  and  supported  by  complimentary  experiments 
in  particmar  cases. 
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The  following  types  of  aircraft  are  under  consideration  or  may 
be  of  interest: 

1.  Flying  platforms  or  "jeeps". 

2.  Lifting  propellers  in  wing  or  fuselage. 

3.  Direct  jet-lift  and  propulsion. 

4.  Combination  of  lifting  propellers  and  jet  propulsion 
with  means  for  conversion  during  transition. 

The  aerodynamic  problem  which  is  the  restricted  subject  of  the 
present  investigation  may  be  divided  into  the  following  subjects: 

1.  Propeller  or  fan. 

2.  Duct  design. 

3.  Problems  of  combination  of  propeller  and  duct. 

4.  Baffles  and  inlet  vanes. 

5.  Interaction  of  ducted  propeller  and  wing  or  fuselage. 

6.  The  effect  of  the  jet  on  the  lower  surface  of  wing  or  fuselage. 
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It  is  conceivable  that  one  or  more  of  the  four  types  of  aircraft 

indicated  above  may  arrive  at  a  practical  stage  and  become  generally 

operational.  In  the  meantime,  it  is  essential  to  arrive  at  the  basic 

design  philosophy  to  the  extent  that  ail  the  fundamental  aspects  of  the 

problems  are  understood.  The  aerodynamic  considerations  are  covered 

( 

under  the  above  headings  as  other  problems  are  of  a  nature  common  to 
any  aircraft.  We  shall  briefly  cover  the  listed  subjects  and  give  reference 
to  the  appropriate  treatment  given  in  this  presentation. 

1.  Propeller  or  Fan 

The  propeller  design  problem  represents  a  rather  straightforward 
case,  except  that  the  propeller  in  a  duct  must  comply  with  the  design 
instruction  given  by  the  formulas  in  Section  B  of  the  present  report.  In 
contrast  to  the  case  of  the  normal  type  propeller  employing  the  Goldstein 
load-distribution  or  the  contra-rotating  type  with  the  distribution  given  by 
the  author,  the  prese  t  case  of  a  propeller  in  a  circular  duct  requires  a 
radial  loading  which  increases  towards  the  tip.  Tables  for  this  case  are 
not  yet  available.  However,  anyone  skilled  in  the  art  of  the  propeller 
theory  can  obtain  approximate  graphs  for  the  loading  of  multiblade 
propellers  in  a  circular  duct  based  on  the  formulas  of  Section  B. 

Design  Method:  Comply  with  the  condition  of  ideal  disk  loading  given  in 
Section  B.  (Tables  not  yet  available.  ) 
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2.  Duct  Design 

Duct  streamlines  for  an  ideal  duct  are  given  in  Ref.Al,  Volume  I, 

Chapter  IV,  and  in  Figures  I  and  II,  Pages  72  and  73.  It  may  be  possible 
to  design  an  ideal  duct  in  some  cases.  The  ideal  duct  is  defined  as  a  duct 
whose  inlet  wall  conforms  to  any  one  of  the  streamlines  shown  in  these 
figures.  Such  a  duct  as  distinguished  from  any  other  duct  provides  a 
constant  axial  component  of  the  inlet  velocity  at  the  plane  of  the  propeller. 

Any  other  duct  design,  for  instance  the  bellmouth  duct,  causes  deviations 
from  the  ideal  case  by  producing  an  excess  velocity  at  the  circumference 
or  at  the  tip  of  the  propeller. 

Design  Method:  Employ  the  ideal  streamline  shape  when  possible  or 
obtain  radial  velocity  distribution  in  any  other  case  by  calculation  or,  if 
necessary,  by  a  model  test. 

3.  Problems  of  Combination  of  Propeller  and  Duct 

The  ideal  cases  of  propeller  and  duct  defined  above  fit  together 
perfectly.  If  deviations  are  necessary,  as  is  usually  the  case  in  the  duct 
design,  calculate  the  resulting  velocity  distribution,  which  normally 
exhibits  an  increase  in  the  velocity  near  the  circumference.  Such  is  the 
case  both  for  short  ducts  and  for  bellmouth  ducts.  This  work  can  only  be 
done  by  engineers  skilled  in  the  art  and  has  to  be  done  as  each  case  will 
be  different.  Experimental  models  or  the  electric  analogy  method  may 

also  be  employed.  In  designing  the  propeller,  make  the  twist  distribution  |j 
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comply  exactly  with  the  particular  axial  velocity  distribution,  but 
maintain  the  loading  as  in  the  ideal  case. 

Design  Method:  Match  propeller  with  particular  imposed  axial  velocity 
distribution  of  the  duct  proper. 

4.  Baffles  and  Inlet  Vanes 

As  the  propeller,  under  no  circumstances,  will  operate  efficiently 
with  a  non-uniform  or  skewed  flow  distribution,  there  are  obviously  cases 
in  which  baffles  or  vanes  of  some  type  should  be  employed.  While  outlet 
vanes  may  be  used  for  control,  no  particular  problem  of  efficiency  is 
involved.  In  the  design  of  inlet  vanes,  it  is  imperative  to  employ  the 
existing  potential  theory.  (See  reference  in  present  volume,  Section  G.  ) 
Tests  performed  as  part  of  the  present  study  confirm  theory  and  show 
remarkably  high  efficiency.  (Sections  H  and  I.  ) 

Design  Me ■  rod:  Design  baffles  according  to  potential  theory  as  indicated 
in  Section  G.  In  complex  cases  experimental  work  may  be  needed. 

5.  Interaction  of  Ducted  Propeller  and  Wing  or  Fuselage 

This  aspect  is  of  importance  in  regard  to  the  required  control 
forces,  lift,  drag,  and  general  performance  of  the  aircraft.  The  first 
item  of  concern  is  the  pitch-up  moment.  Numerical  values  have  been 
given  in  Ref.Al,  Volume  I,  Chapter  II  for  a  sink  in  a  rectangular  plate 
and  in  Chapter  III  for  the  classical  case  of  a  two-dimensional  line  sink 
in  a  wing. 
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Design  Method:  To  obtain  pitch-up  moment,  employ  the  simple  formula 
given  in  Ref.  Chapter  II.  The  pitch-up  moment  arm  is  given  in 
Chapter  III  under  the  formula  III— Z 1  • 

There  is  also  a  lift  increase  due  to  induced  circulation  on  the 
wing  which  can  be  estimated  with  sufficient  accuracy  by  formula  III- 18 
of  the  same  chapter. 

6.  The  Effect  of  the  Jet  on  the  Lower  Surface  of  the  Wing  or  Fuselage 

This  subject  is  of  more  empirical  nature  and  represents  the  only 
area  which  cannot,  for  the  most  part,  be  covered  by  theory.  The  only 
exception  is  the  direct  displacement  effect  of  the  jet  on  the  pressure 
distribution  for  the  lower  surface  of  the  wing  or  fuselage.  This  subject 
is  treated  in  the  present  paper  in  Section  E.  Due  essentially  to  the  high 
velocity  on  the  lateral  sides  of  the  jet  there  results  a  negative  lift  in 
magnitude  equal  to  the  product  of  the  velocity  head  times  the  area  of  the 
jet  cross  section  and  a  negative  lift  coefficient  approaching  the  value  of 
unity,  if  the  adjacent  flat  area  is  large  enough. 

Design  Method:  Estimate  extent  of  adversely  affected  area  and  employ 
negative  lift  and  thrust  coefficients  given  in  Section  E  and  shown  in 
Figures  El  andE2. 

There  is,  however,  a  much  larger  negative  contribution  both  to 
the  lift  and  to  the  adverse  pitching  moment.  This  subject  has  been  treated 
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in  the  present  volume,  Section  F.  The  deficiency  in  pressure  recovery 
which  is  the  cause  of  the  drag  of  a  cylinder  and  which  exists  behind  the 
jet  produces  a  pressure-deficiency  on  the  lower  wing  surface  behind  the 
jet.  Tests  analyzed  in  Section  F  indicate  that  there  exists  a 

drastic  decrease  of  lift  and  a  large  increase  in  the  pitching  moment. 
However,  since  the  drag  of  a  cylindrical  body  is  reduced  to  a  fraction 
of  its  value  beyond  a  certain  critical  Reynolds  number,  it  is  probable 
that  the  effect  is  less  serious  than  indicated. 

Design  Method:  Yet  inadequate.  Tests  of  pressure  distribution  at 
higher  Reynolds  number  (R>  10  )  are  needed.  In  the  meantime, 
employ  values  of  one  quarter  of  the  ones  indicated  from  small  scale 
tests  in  Section  F. 

There  exists  a  few  cases  in  which  the  theoretical  treatment  is 
either  too  laborious  or  inadequate.  In  conjunction  with  a  design  of  any 
consequence,  certain  complimentary  experiments  could  be  carried  on  to 
considerable  advantage. 

The  first  problem  concerns  the  matching  of  the  propeller  to  the 
duct.  A  simple  experiment  may  be  performed,  the  purpose  of  which  is 
to  establish  the  flow  pattern  of  the  duct-inlet  at  zero  forward  velocity. 
This  test  is  most  conveniently  performed  by  attaching  a  cylindrical 
extension  of  some  length  to  the  model  duct  outlet  to  remove  the  local 
effect  of  the  suction  device.  Measure  the  velocity  distribution  along  the 
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radius  of  the  duct  at  the  proposed  plane  of  the  propeller.  This  is  also 
an  excellent  case  in  which  the  electric  analogy  may  be  employed.  The 
advantage  is  that  the  effect  of  the  viscosity  or  Reynolds  number  is  then 
completely  eliminated.  Thus  a  very  small  model  may  be  employed  to 
represent  properly  the  condition  of  the  actual  flow  at  a  high  Reynolds 
number.  The  matching  propeller  can  then  be  designed.  Testing  of  a 
model  propeller  will  not  be  necessary.  The  model  duct  should  be  of  at 
least  one-foot  diameter  to  insure  reasonable  Reynolds  numbers  in  a 
typical  low  speed  wind  tunnel. 

The  second  experimental  problem  concerns  the  empirical  effect 
of  the  obstruction  of  the  jet  on  the  flow  pattern  on  the  lower  side  of  a 
wing.  As  has  been  pointed  out,  the  deficiency  in  pressure  recovery 
behind  the  jet  causes  a  loss  in  lift  and  creates  an  additional  unfavorable 
pitch-up  moment  of  some  magnitude.  In  addition,  the  jet  exhibits  a 
pumping  action  on  the  lower  side  of  the  wing  tending  to  aggrevate  these 
effects.  To  obtain  numerical  values  of  the  loss  in  lift  and  the  adverse 
pitching  moment  and  to  investigate  means  for  improvements,  a  fan-in¬ 
wing  model  test  should  be  run  at  a  reasonable  Reynolds  number  and  the 
model  should  be  equipped  with  means  for  measuring  the  pressure 
distribution  on  the  lower  surface  as  affected  by  the  jet  in  the  range  of 
forward  speeds  corresponding  to  transition. 
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SECTION  B 


THE  IDEAL  LOADING  OF  SINGLE-ROTATION  PROPELLERS 
IN  A  CIRCULAR  DUCT 


I.  INTRODUCTION 


In  this  problem,  the  flow  field  of  a  single  shrouded  propeller  is 
examined.  The  shroud  is  considered  as  an  infinite  cylindrical  surface, 
and  the  velocity  potential  function  (j)  which  describes  the  flow  within 
the  shroud  is  determined.  Use  has  been  made  of  Sydney  Goldstein's  method 
in  solving  the  boundary  problem  for  the  ideal  case  of  a  single  rotating 
propeller. 

In  this  case  a  two-bladed  propeller  will  be  considered  and  the  helical 
surface  produced  by  the  propeller  is  given  by 


0  -  - =  0  c/u  9r 


CD 


where  r,  0  and  z  are  the  cylindrical  polar  coordinates  with  the  helix 
axis  as  the  reference  axis.  The  axial  displacement  of  the  spiral  is  w 
and  the  expression  for  the  slope  at  any  radius  r  is 


'tosn  c< 


__  V-f-  W 

n-d®  s\,  go 


(2) 


Thus,  the  velocity  normal  to  the  surface  is 

9  cj> 

w  CCHD  CK  =  — —  co<L-  °< - — r  (3) 

o  ^  /l9  y 

Tne  above  boundary  condition  is  given  for  the  helix  surface  only  for  r  <  R 
where  R  is  the  radius  of  the  helix.  The  boundary  condition  for  r  a  R 
is  that  3  <j>  /3r  “  0. 


9 


Changing  the  variables  r,  0  and  z  into  p,  and  J  given  by 

and  taking  into  consideration  that  J  °  constant;  gives  a  spiral  line, 
the  three-dimensional  problem  reduces  to  a  two-dimensional  one.  From 
Eqs.  (It),  one  obtains  the  relations 

li.  -  _  1:0  a  4>  gi  gi. 

3^  ”  v+w  3  ^  ae  "  5 

Making  use  of  Fq.  (3),  we  find  the  following  for  the  boundary  conditions  in 
the  (p,  ^  )  system 

3(f)  ^  V+w  ,  a<l>  .  +  R 

3  $  i-^2-  CO  T  3/*  /  ^  0  V+W  w 

Substituting  the  following  function  for  the  velocity  potential 


=  (j) 


the  boundary  conditions  become 

3^  A  25Ll  -  n  ,fi1 

3c)  1+/- 

The  differential  eouation  for  the  potential  in  conventional  cylindrical 
coordinates  is  given  by 


i  a  /  9?  \  ,  3Z^  .  ll£_  _  n 

v  ^  -  /<-  r  sa)+a,2562  3^  0 


From  Eqs.  (h)  and  *5),  we  obtain  the  following  relations  among  the  partial 
derivatives  of  the  function  (j>  in  the  two  coordinate  systems: 


36  =  3^ 
32  _  31 

ae2  = 


9  _____  co  9 

3 1  V+w  d  ^ 

3Z  =  /  co  Y 
\v+wy  a^2- 
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3a  V-hw  3/^ 


Substituting  for  the  above  relations  in  Eq.  (9),  we  find  the  Laplacian  in 
the  transformed  (p,  ^  )  coordinate  system 

fg^)  ■*- 0 +/~‘)  fp  0  o 

Thus,  the  problem  of  the  flow  field  of  the  shrouded  propeller  is  given 
by  the  solution  of  the  partial  differential  Eo.  (10)  with  the  boundary 
conditions  given  by  Eos.  (8). 

II.  SOLUTION  OF  THE  POTENTIAL  PROBLEM 

Since  it  has  been  assumed  that  the  shroud  is  a  circular  cylindrical 
surface  extending  to  infinity,  we  are  interested  in  obtaining  only  the 
flow  field  in  the  region  r  <  R.  Changing  the  dependent  variable  into 


= 


the  above  Ec .  (10)  becomes 


(' =  p;)  (fg’)  S  08 

From  Eos.  (8),  the  boundary  conditions  for  the  function  are  given  by 

^---0  4  ^-1  =  0  (13) 

Expanding  the  independent  variable  ^  on  the  "right-hand  side"  of  Eq.  (12)  in  a 
Fourier-series  for  the  interval  0  «<  ^  •*<  it  we  find 

\  =  x  -  ±  i  (V) 


Fourier-series  for  the  interval  0 


<  it  we  find 


tj.  CoAs  (2  +  i)  ^ 

F  £=0  7^0* 


II 


Assuming  a  solution  ^  which  is  expanded  in  the  form 

%  =  -jo  O)  +  jt  W  «*-  (*-»  + 1)  5  (W) 

and  substituting  both  Eqs.  (Hi)  and  (l£)  into  Eq.  (12),  we  find  the  following 
by  equating  the  coefficients  of  cos(2m+l)  \  > 


Since  the  potential  ^  is  finite  at  r  ■  0,  the  above  Eq.  (16)  gives 


Using  for  the  function  f  (p),  the  expression 

m 


J  /»=__£ _ 1 _ 

1  7r  (2^+lf  ) 


1+/*- 


O) 


and  substituting  into  Eq.  (17),  we  find  the  following  differential  equation 

for  g  (p). 
m 

/  ^  \  2  ^  ^ 


A  particular  solution  of  Eq.  (20)  is 


X 

L(?)+ye 


where  t  (z)  is  given  by  the  series 

1)V  ,  j 


71^  "i  ^  * 


2*V  (?~V%W 


+  -•  + 


(zW)  yz) 


r-rrf-(  22) 


and'  v  “  2m  1-1  &  z  =  (2m-t-l)p. 
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We  define  a  new  function  T  ( [2m+ 1\  p. ) . 

l,2ra-t-l 


%SMr>  (23) 

Then  the  above  function  T  is  also  a  particular  solution  of  Eq.  (20). 

1,2m +  1 

The  general  solution  of  Eq.  (20)  which  remains  finite  on  the  axis,  (i.e., 
p.  “  0)  is  thus  given  by 

(/*)=  + b~  w 

Substituting  into  Eq.  (11)  for  from  Eq.  (lii),  for  Cp  from  Eo.  (15), 

for  f  from  Eq.  (18),  and  for  f  from  Eq.  (19),  we  find 
o  m 

Tidbit 


-  r  Ji  +  a 

OW'J’r  (z-w.fi)2-  ~  T  ([2^-mI*. 


Since  the  function  (j>  is  an  odd  function,  it  is  obvious  that  for  the 
interval  -tt  <.  ^  0,  we  find  the  same  expression  for  (f  with  the 

opposite  sign.  At  ^  -  0  or  tr,  the  potential  function  ^  has  a 
discontinuity  which  is  given  by 


(2a*+i)5 


r^i  g  4  Jz 

L?V  7r  £o  (z^,y  ^ 


The  unknown  coefficients  a  are  determined  from  the  second  of  the 

m 

boundary  conditions  given  by  Eq.  (8).  Differentiating  the  above  Eq.  (25) 


with  respect  to  p.,  we  get 

CO  f  ,  I  7 

__  I  jf__  '**  MZ^r I 

nr  (2/wti)1  I"  ft 

'  ■  2.^+1 


i/ 

i-  MHH  (27) 
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Making  use  of  Eq.  (27),  the  boundary  condition  given  by  the  second  of  Eqs. 
(8)  gives 

(i^H)z  I'  (2 

Z  +  j 

Substituting  Eq.  (28)  into  Eq.  (25),  we  get 


r  4  * 

9s  * 


i  |T  /r  \  t;L+i(Cwh:W 


l!  ([2/w+iLmJ  Z^+l 

ZA«-h|  x 


Thus,  the  velocity  potential  <j)  from  Eq.  (7)  becomes 

4>=  £  jr  ((^) 

VF  <XT  rtr\  ZO  I  l}Z'y*  +  l  '  I  ([Z^fflcy  Z«+|  ' 

1  2.a*+|  ' 


29) 


co^ 


(z^-rl)‘ 


(30) 


Similarly,  the  jump  of  the  velocity  potential  at  |  *  0  or  it,  given  by 
Eq.  (26)  becomes 


&]« 


gv/(Yw)  £  ]  ) 


TCCO 


(31) 


/w^r  0 


III.  VELOCITY  FIELD 


The  velocity  field  within  the  shroud  is  formed  by  the  gradient  of  the 
velocity  potential,  i.e.,  v  ■=  v  (j)  •  Making  use  of  Eqs.  (ii)  and  (5),  we 
find  the  following  for  the  velocity  components  in  the  cylindrical  coordinate 
system 


9<j>  9<j) 

U6=  ZdQ  >  “l  ~  (32) 

In  the  (p,  ^  )  system,  the  above  equations  become 

co  d$  _  co  I  3<f>  _  _co_ 

CA~  v+  w  9/*  }  Vtw  ^  ~3^~ }  §~  Vi-w  3  ^  ^ 


di> 

(A  =  ~ 
/L  3/h 
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From  Eos.  (33),  it  is  obvious  that 


=  V*  “e 

Differentiating  Eq.  (30)  with  respect  to  ^  and  p,,  we  get  respectively 

3<f>  _  4vj(j+w)  •£-  \Tfc  il  A  T  fr>  n 

jco  £■„  (2£,  W  j/  mfl]/>\~(z:+ly  ok ) 


9(j)  4  wfVVw)  ®  /  T  (C2^+i]m^) 

~  21  iT  (Lz^tllu)-  - -r-I  ([2WjV<) 

l  Z^ti  ' 


C0<L(2/vn  tl)^ 

~0-+i)  (35) 


From  Ers.  (32),  (3ii )  and  (35),  the  velocity  components  u  ,  u  and  u 

**  0  2 

become 


w"fc  - 


<rr,  r  D  **  ^  2/WI  t  t 


J  00 

W'  = - 2  A  (>W 

*  ^  ^To  - A 


^2.  /wi  f  ^ 
2  ^  +  / 


V  Z  *>.y> 


where: 


A  (/\  y^)  -  I  (&™i-  /]/*)  — kk?..±j — _ -  j  (jj2./rt\  +  |]u^ 


Al0.y)=  T'  (&-...»-  J&-TT  I'  a-M 

I  ([Z/W-HjMj  Z/ywtl  '' 

2/wt-M  / 
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IV.  LOADING  FUNCTION' 


The  loading  function  K(r)  is  a  dimensionless  quantity  defined  by 

K(a)  r  “XLH  (111) 

where  p  is  the  number  of  blados  or  helix  surfaces  and  f“  (r)  is  the 
circulation  which  is  equal  to  the  potential  difference  across  the  helix 
surface  at  the  radius  r. 

From  Eq.  (31),  the  circulation  f  (r)  becomes 

2_  '  v  (t  j  ft*) 

^  Jt  CO  ^3  0  O'"*1 )  ll^-t  I  X  '  \ 

1**  tl  J 

Introducing  the  dimensionless  quantity 

R  ft3> 

and  making  use  of  Eq.  (Ijl),  we  get  the  following  from  the  above  Eq.  (i|2) 
a  “if  T/(^'^+i1a<)  1 

Km  =  ~  Z  r-~r  t  fs&  f-T  (c^+3.)  (u.) 

y  --2-  (2or,  +  |)  +  ,  j  0>w+(UO  2/*'+l  \ 

b  2^tl  J 

where 


"  V-t 


60  03  R 

/L*  —  —————  £ 

V  +  \/+w 


The  momentum  within  a  volume  confined  between  two  successive  vortex 


surfaces  become 


where  f~  is  the  circulation  function  for  one  blade  and  A  is  the  area  of 
the  projected  cross-section  of  one  turn  of  the  spiral,  which  is  equal  to  the 
cross-section  of  the  shroud.  Introducing  again  the  dimensionless  quantity 
't  ,  the  momentum  per  surface  per  turn  becomes 
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fiTri 

M=  ^w^VR2  ± 

/ o 

Tlius,  the  mass  coefficient  for  single -rotation  propellers  is 


K(T>  fytoitcLQ- 


(1*7) 


w  - 


7r 


27r 


I 

K(?> 6)toL'toL  0 


Kfr)  TdT 


(1*8) 


Jo  ->  a  */o 

Substituting  for  K('T)  from  Eq.  (hi*)  into  the  above  equation,  the  mass 
coefficient  becomes 


j 


(1*9) 


From  Eos.  (hi*)  and  (1*9)  for  the  loading  function  K  and  the  mass 


coefficient,  ^  ,  we  find  the  design  relations  of  a  propeller. 
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V.  APPENDIX 


So  long  as  p,  is  not  'boo  small,  we  might  make  the  following  approxima¬ 


tions.  The  functions  T  given  by  Eq.  (23)  are  approximated  by 

l,2m-f-l 


T 


l}Zsrr\  +  l  /J/* 

Thus,  its  derivative  will  be 


(A-l) 


T 


/ 


'  I  *-h)  Zsr>\-hl  (jptA  ) 
From  the  approximate  relations 


(2^-t-i  )X1/vv,+j  _  K 

- - - -  Z  sr\  -~ 

I  K 

2.  rv*\  f  | 


I*  . .  .  / 1  .  I  ^ 
Zs* 


■z  v/  + 


Z/W\  •  I  -f-  2/n 


’) 


l< 


^  -  /itX 


we  find 


I' 


1 


^  y/t/* 


(A-2) 


(A-3) 


(A-l.) 


Solving  the  above  eauations  for  I  .  we  find 

f  2m -f-1 

I  .  ^  e 

2/v^ +j 


and  its  derivative  becomes 


I 


t  *  e 

Z^tl  ' 


■  j/*  z  ■+  A 


(A-*) 


(A-6) 
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Equations  (A— 5 )  and  (A— 6)  are  approximated  by 


r  *  e 

2'vn  +  | 


2wtl 


(A— 7 ) 


Zm)f! 


I '  *  •Jiiy?  e 

E'm-H  ' 


(A-8) 


Making  use  of  Eqs.  (A-l),  (A-2),  (A-7)  and  (A-8),  we  find  the  following  for 

the  functions  A  (p.  ,|i)  and  A1  (p  ,p)  given  by  Eqs.  (39)  and  (1*0). 
mo  mo 

A  />  - ! - f££_  e _  (A.o) 


A,  (a/)*  ^ 


■w  ^  I  I  ■  Ml—!  Pi 

TV* //>“.*■ 


(A- 10) 


Substituting  Eqs.  (A-9)  and  (A-10)  into  Eqs.  (1*1*)  and  (1*9),  we  find  an 
approximate  estimate  for  the  loading  function  K  and  the  mass  coefficient 
nj.  of  ti*C  propeller. 

Similarly,  Eqs.  (36),  (37)  and  (38)  will  give  an  approximation  for  the 
flow  field  within  the  shroud. 
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SECTION  C 


FLOW  FIELD  FOR  NON-UNIFORM  SINK  DISTRIBUTION 
IN  A  CIRCULAR  DISK 

I.  INTRODUCTION 

In  the  Interim  Report  Number  7  of  contract  DA  I4J4  —17 7 — TC  — 606 ,  the  flow 
field  for  uniform  sink  distribution  in  a  circular  disk  was  derived.  In  this 
paper  the  case  of  non-uniform  sink  distribution  will  be  examined.  For 
simplification,  axially  symmetric  cases  will  be  considered  and  a  step-wise 
constant  distribution  over  circular  rings  will  be  assumed. 


The  cas  cf  a  uniform  sink  distribution  in  an  annulus  is  found  as  a 
special  case,  by  considering  the  strength  of  the  sinks  in  the  central  region 
as  zero. 

The  flow  field  of  a  propeller  in  an  annulus  of  an  infinite  plate  is 
then  represented  by  a  uniform  sink  distribution  at  the  annulus  and  a 
superposition  of  a  uniform  flow  in  the  serai-infinite  region  of  the  cylindrical 
annulus  behind  the  propeller. 

Finally,  the  case  of  a  continuous  and  axially  symmetric  sink  distribution 
is  examined  which  is  expressed  as  a  polynomial  of  even  powers  of  the  radius  rQ. 
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II.  VELOCITY  POTENTIAL 


For  the  case  of  only  one  step  in  the  distribution  of  the  sinks, 
we  have  the  following 


Figure  C2 

From  the  velocity  potential  of  a  point  sink  of  strength  <f  Q  and 
an  angle  co  »  90°  at  a  distance  rQ  from  the  origin,  we  have 
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In  Figure C2 the  step-vise  sink  distribution  can  be  considered  as 
the  algebraic  summation  of  two  uniform  sink  distributions  over  two  co¬ 
axial  disks  of  different  radius 


Thus,  in  the  above  case,  the  velocity  potential  is  found  as  follows: 


Region  A :  r  ;> 


For  r>  Rp  i.e.,  r/rQ  >  1,  the  velocity  potential, for  the  sink 
distribution  of  strength  Qp  along  the  axis  of  symmetry  0  ■  0  is 
found  by  integrating  Eqs.  (l)  from  0  to  R^,  i.e., 

<£,  M  =  f  '  J  Z.  Pr  (a)^drc  d(t 

6~0)cj~cIo>  A  A  7  n=c 

ft 

Integrating  and  substituting  S  x  tt  R^£  ■  Q^,  we  find 

h  fro W2l|  JStl 


(2) 


(3) 


wnere 


Q,  =  zrRfA 


(h) 


Similarly,  Tor  the  second  uniform  distribution  cfQj  from  0  to  RQ, 


we  gel 


L  (no)  - 


'o  /o 


c/~Qz 


1(^-Jpn(o)r0^^ 

r\~o  '  1 


(5) 


Integrating  and  substituting  CL  a  cT  Q  x  tt  R  A  we  find 

C  4  O 


x  r  \  02_  y-  ft  (°)  ^ph-_ 


(6) 
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where 


Q2  =  27rR*\  =  zrf?/(v,-V0) 


Thus,  for  r  >  IL  the  velocity  potential  along  the  axis  of  symmetry 


becomes 


^  £  4^-  Or (t^J 


Since  the  flow  is  axial  symmetric,  the  general  solution  of  the 
Laplace  equation  is  given  by 

$  (r,  e)  =  £_  Anrnpn  (C~e)+l  eK  -4+t  P„  (co*fl)  , 

n:o  r 

It  is  obvious  that  the  velocity  potential  in  the  region  r  >*L  f0r  any 


angle  0  is  given  by 


i  (p,e)=4r£ 

V\v>  27TV  TTo  >°  +  2-  ^ 


^  (<^e) 


tor  r  <  r0  <  Rc 


Region  B :  r  <  R0 

For  r  <  R„  we  have  that 
^  o 

ir/r  >  1  for  o  <  ir ■  <  r 

'  (11 
r/r0  <1  r  <  Tc  <  f?0 

Thus,  the  velocity  potential,  for  the  sink  distribution  of  strength  £  Q^, 
along  the  axis  of  symmetry  9  ”  0  is  found  by  integrating  Eqs.  (1)  from 
0  to  r  and  from  r  to  Rj,  respectively,  i.e.. 


*  Y'  ■'0 
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Performing  the  integration,  we  find 


1  '  '  iTrRf  n=o  h-t  n+z  l  Ru 


where 


40  co 

Z 

n=o 


Similarly,  for  the  second  sink  distribution  of  strength  <TQ 
we  get 

-  — 

*2  '  ’  ^  27rF?0a  h=o  n"1  ntt.  \  R«  / 

Thus,  for  r  <  Rq  .  the  velocity  potential  along  the  axis  of  symmetry 


becomes 


i  /  \  r  r  1  P»  (o)  Q<  r  2h-f1  /  r  Y1"* 

^g(1'.0)- $|  27f  Z  h„1  j  R(i  nti  (  rJ  J 


2  Htl 

fV 

htZ 

-(£)" 


X  (r  _  c:*  &  R)  C°)  )  2r\Jr\  Q 1  Qz  L  Q2.  /  r  Y1”^  ^  /  r  Y"1"^  1 

V  ,w  -  ar  ~  VSTT  ‘ '  [W  "  W  J 


where 


Qj  _  Qz.  _  zTrR 


r,£  r? 


R,z  ’< 


mM  -  27TV0 


Kj  f'o  1^1  1 'O 

Making  use  of  Ea„  (9),  since  the  flow  field  is  axially  symmetric,  we  get 

A  -  _L  (Sl.  ®iJ\  •? *  ftW  Hii 

/M  ~  OTr\  0.2  P.w4.  U-»H  K.+.Z. 


k-1  K-ci 


1~  Z7T  V,  R,2  f?oV  k^o  k'1 

_L  2&(Ji L_  _  _2sJ\ 

h~  27r  n-1  V  R.nt{  K.”*y 


A _ i_  gw  Q 

“h~~.29r  n-1  \(?.' 


6  50 


25 


(19) 


From  the  sum  (see  Interim  Report  No.  7) 

jr*  2k+1  _  _j 

K*o  K-1  K« 


and  from  Eqs.  (9)  and  (18)  we  have  that  the  velocity  potential  in  the 
region  r  <  RQ  for  any  angle  0  is  given  by 


PnO>) 

(n-o 


(20) 


Region  C :  RQ<  r  < 


For  RQ<  r  <  R-j_  we  have 

that 

r/r0  >  1 

for 

R0«C  r0<  r 

r/rQ  <  1 

for 

r<r0<Ri 

Thus,  the  velocity  potential^for  the  sink  distribution  cTQ-^,  along  the  axis 
of  symmetry  G  ■  0,  is  found  by  integrating  Eqs.  (1)  from  0  to  r  and 
from  r  to  R^,  i.e., 


Performing  the  integration,  we  get 


Similarly,  for  the  second  sink  distribution  cTQj,  we  have 


i  (r, o)  = 


2  7T 


cTQ; 

HTr 


r  p^o)  r°  dr°  dcf 


(23) 


or 


I  /  \  |  f?(°)  R0n 


*=o  r""’  (2ii) 

Thus,  for  Rq<  r  <  the  velocity  potential  along  the  axis  of  symmetry 
becomes 


<r^r  P  (°) 

11  " 


2n»1 

n-t-L 


Qz  £  R<£)  R," 
frt  n+z  r"+< 


(25) 


or 


\  (r>°)  = 


Q) 

Z'Trf?.] 


_r  \  R,  (°)  2  n-H 

^*1  '  *zO  h_1  h  +  2- 


+- 


C  26) 


Making  use  of  Eq.  (19),  the  above  equation  becomes 


Qt 

/JL. 

\ 

^7rR| 

J 

27rR\ 

£* 
n  «o 

r  £(o) 
n-1 

(27) 

Qz 

co 

n-o 

P  (»)  ( 

'(?„  V11"1 

27rR0 

—J 

27 


I 


Thus,  the  velocity  potential  in  the  region  RQ<  r<R^  and  any  angle  0 
becomes 


%  ^(o)  /  r  \n  Qz 

2-o  n-1  (r,.)+29tR0 
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n-o 


50) 

n  4-2 


III.  VELOCITIES 


(28) 


The  velocity  components  for  the  three  regions,  A,  B  and  C,  are  given 


from  the  derivatives 


u 


r 


dr 


U 


(29) 


Region  A:  r>Rl 


From  Eos.  (10)  and  (29)  we  find 


00 


ut~~  Zirr.*  h+k 


=  -  J--1  ^  - Qt  (-T-) 


<^6“) 

ot/U. 


(30) 


Region  B :  r  <  R0 


From  Eqs.  (20)  and  (29),  we  find 
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Region  C :  RQ  <  r  < 


From  Eqs.  (28)  and  (29),  we  find 


U  =- Sk-h.  CO<u 0  - 


2]_^*-ngW/jr_nj2t_S  HLpfjSsS*2 

nTo  n"'(  'RW  27rP0X  £~o  h+l  "  ^r/ 


WG  I  Ql  [Ql  S^/r  vh  ii.  I /.Sir1]  *2^ 
e  2^  I  R,2  t(_R|a-  «.  n-1  R02  „4  »U  (.  r  )  J  0^, 


where 


0^  -  27T  f?|  Vj 

Qz--  2*f?0s ?(vf-V0) 

Si-  -  -§r-  -  2?TV 
R7  R*  fV° 


£  (<**©) 

(32) 


(33) 


IV.  STREAMFUNCTION 


From  the  continuity  ecruation  in  spherical  coordinates,  i.e.. 


1 


ar  (*-l“r)+F3=9  7i r(“e^«)=° 


the  streamfunction  is  defined  as  follows 

at 


ar 


=.  -  2.9T r  Uq  0^0 


=  2Vrrzur 

So  that  the  flux  is  given  by 

at  =  + 

■  ar  30 


(3U) 


(3« 


(36) 
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Region  A:  r>R1 


From  Eqs,  (32)  and  (30)  we  have 

|i-  =  -  iTt  we  A-e  = 


1  2 
=  —  /w*  £ 


Integrating  we  get 


Differentiating  the  above  equation  with  respect  to  9,  we  find 


2  Qc^a  ^  .^q 

d-^  c^u 


From  the  second  of  Eos,  (32)  and  from  the  first  of  Eqs.  (30),  we  have  that 
-|~  =  27Trz  urM~Q 

r  (W 

=  ->“&£  P„  (°)  Qi  (t1)  -  P2  (t-)  £  f/0 

Comparing  Eos,  (i|0)  and  (hi)  and  rearranging  terms,  we  get 


Since  the  differential  ecruation  in  the  brackets  is  equal  to  0  the  above 
equation  becomes 


30 


/'(»>" 


/O^  0 


Integrating,  we  get 


-/ (e)  =  (<** .e+q)  , 

where  the  constant  is  evaluated  from  the  oondition  +  (0)  -  0  for 


0  ■  0.  Consequently, 


c,=  -< 


And  Eq.  (38)  for  the  streamfunction  in  the  region  r  ?  R^  becomes 


Region  B :  r  <  R 


We  follow  the  same  procedure  as  that  for  the  Region  A, 


From  Eqs.  (35)  and  (31),  we  have 
«  ~  -  —  Z.7T 0 


=  -rA^ 


-aJ/Qi  °2  V  <?*■  ^  ^  3Lf  h  \n_1 

'p‘(rTJ 


Integrating,  we  find 


f  =:  -rVl 
'B 


Differentiating  the  above  eauation  with  respect  to  0  we  get 


,  _Z  •  3 
-j-  r  xuv- 


0=0  (nH)(htl)  olu2  0 
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From  the  second  of  Eqs.  (35)  and  from  the  first  of  Eqs.  (31),  we  have  that 


Comparing  Eqs.  (1*9)  and  (50)  and  rearranging  terms,  we  find 


r*  ■  a  4r*  pn(°l_  \SLfJLn 

r^6h  MM)  Rel\RoJ  \[  a) 


Since  the  differential  equation  in  the  brackets  is  equal  to  0,  the  above 


equation  becomes 


Integrating,  we  get 


/(e) = 

^(e)=c 


where  the  constant  C^,  is  evaluated  from  the  condition  '"f  (0)  ■  0  for 


0=0,  Consequently, 


2(e>c,= 


And  Eq.  (1+G)  for  the  streamfunction  in  the  region  r  <  Rq  becomes 

y(r6)..r Wfl life L - S*-V Zl°2  *5.1 

rB(r,e)-  t-A-e  R^fUJ^) 


Region  C :  RQ  <  r  <  R^ 

From  Eos.  (35)  and  (32)  we  have 
3Nfl 

-—---ZTrr 


—  "T Q 


|q,  .fa,  pn( o)f  r  Qt  »  pm  /R.r1!  oL?k 
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Integrating,  we  get 

»0f<5f  Jo,  g*  RW  |rf  gt  2  RW 

"c  ®|2f?|i+[f?lJ'Sv  Si  n(n«)(  *•/  J  <^u 

^  - 

■hj(S) 


(57) 


Differentiating  the  above  equation  with  respect  to  0  we  find 


2/&/-*-0c/w.0  + 


(58) 


+rz<o*^6 


MMlRj  R 


n-l  Qj_ 


z 


P»  /jvr 2 

h(h-f^)  \  r  / 


or 


9± 

ae 


2  %C<0  /Ro\nn1 
r  /  J 


/oI^/@  c<^0  +^;(q) 


(59) 


From  the  second  of  Eos.  (35)  and  from  the  first  of  Eqs.  (32),  we  have  that 


It 

39 


E7r rz  Up  a3aQ 


(60) 


=  -r2^G 


Co</Q  -r 


2l 
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Comparing  Eqs.  (59)  and  (60)  and  rearranging  terms,  we  get 


Since  the  differential  eouation  in  the  brackets  is  equal  to  0  the  above 
equation  becomes 

■j  (e)  =  0  (62) 

Integrating  with  respect  to  0  we  get 

1-(6)=C3  (63) 

where  the  constant  is  evaluated  from  the  condition  ^  (0)  for  0  *  0* 
Consequently, 


^.(o)=C3=0 


(6k) 


And.  E<j.  (57)  for  the  streamf unction  in  the  region  RQ<  r  becomes 


vu,  nw-r^-r  Jh.  ■?»  p"  fa  f  r.f  Qz.  ®  pn(°)  ,  2L 


J 


y  (65) 


V.  SINK  DISTRIBUTION  IN  AN  ANNULUS 

For  a  uniform  sink  distribution  in  an  annulus  with  internal  radius 
Rq  and  external  radius  R^  the  flow  field  is  derived  by  putting  VQ  a  0 
in  Figure  C£. 


34 


(66) 
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Similarly,  substituting  for  and  from  Eqs.  (66)  into  Eqs.  (1*6), 
(5$)  and  (65),  the  streamf unction  in  the  Regions  A,  B  and  C  becomes: 


Region  A :  r  ;>  R^ 


d 


T 


Region  B :  r  <  R 


Region  C:  R0<  r  <% 


g<R>M 

d^. 


(70) 


(71) 


•r 


+f 


(72) 


VI.  PROPELLER  IN  AN  Am'LUS 

The  flow  field  of  a  propeller  in  an  annulus  of  an  infinite  plate 
might  be  represented  by  a  uniform  sink  distribution  at  the  annulus  with 
strength  eoual  to  twice  the  normal  velocity  at  the  propeller  plane. 

At  the  sink  plane  there  are  symmetric  radial  flow  velocities  and 
a  uniform  axial  velocity* 

Thus,  the  flow  field  in  the  region  above  the  propeller  is  represented 
by  the  Eqs.  (67),  (68),  (69)  and  (70),  (71)  and  (72). 
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In  the  region  under  the  sink  plane,  an  annular  jet  emerges  with  an 
increased  total  energy  due  to  the  energy  added  by  the  loading  of  the 
propeller. 

For  this  region,  a  solution  to  the  potential  flow  field  might  be 
found  by  adding  a  uniform  flow  of  double  the  axial  velocity,  in  the 
cylindrical  annular  region  with  cross-sectional  area  equal  to  that  of  the 
annular  sink. 

Thus,  superimposing  the  annular  sink  flow  field  to  the  uniform 
velocity  field,  we  get  a  flow  field  which  has  an  annular  jet  with  a  contraction 
ratio  equal  to  two. 

At  infinity,  the  mass  flow  from  the  annulus  of  the  plate  is  eaual  to 
the  mass  inflow  from  the  cylindrical  surface  due  to  the  sinks.  A  sketch 
of  the  flow  field  is  given  in  the  following  Fig.  C4. 
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VII.  CONTINUOUS  NON-UNIFORM  SINK  DISTRIBUTION 

For  axial  symmetric  non-uniform  sink  distribution  we  have  that  the 
strength  of  the  sinks  is  given  as  a  function  of  the  radius  rQ,  i.e., 


Figure  C5 


39 


40 


I 


r 

for 

o  <  r  ■<  r 

0 

r 

r 

<  i 

for 

r  <  r  <  R 

o 

(83) 


Thus,  the  velocity  potential,  for  the  sink  distribution  of  strength  q(r0), 
along  the  axis  of  symmetry  9  ■  0  is  found  by  integrating  Eqs.  (1)  from 
0  to  r  and  from  r  to  R_ ,  respectively,  i.e., 


(8U) 


Substituting  for  q(rQ)  from  Eq,  (l$)  and  integrating  with  respect  to  ^  , 

we  get 


'r  N 


$  (r>o)*-z  Z  a  r*  L  Pn  (») 

I  K'o  rt-6  1 

’Ri  (4 


'  zj  ^  pj°> 

Making  use  of  Eq,  (78),  we  find 

$(0°)  =  t(  Z  b„(^-)K  f=Co)«J, 


+  2f  ilz  r°) 


(85) 


(86) 
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Interchanging  summation  and  integration,  we  get 


§(r,o)4Zpn(c)£bKrl  (*T  *&) 

u  K-o  / 


®  M 


i"i£lZ(0>fcb*  I 


(87) 


Performing  the  integration,  we  find 


M 


n  +  K+Z. 


I-W 


n-K-i 


+ 


(88) 


where 


CD  j,  i  ® 

Z*-L 

h=o  n=o 

n  A  K  +•! 


(89) 


For  k  ■  even,  i.e.,  (k*  “  2k),  we  have  an  even  polynomial  for  q(rQ) 
and 


P*.  (o)  -  P  (o)  ~  o 

K>|  W  2  K-H  '  ' 


(90) 


Thus,  Eq.  (88)  becomes 


n’°  K:o 


H+2K+Z. 


+- 


(91) 


And  the  velocity  potential  in  the  region  r  <  R-j_  for  any  angle  0  is 


given  by 

&aM= 


2frR, 


l  wt  7^rA*\4  ^ 


A,  f** 


AL 


27T  Pi  n&O 


^  P.WZ 


(92) 


K*o 
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where 


(93) 


The  velocity  components  for  the  two  regions  are  given  from  Eqs*  (29) 
and  Eqs.  (82)  and  (92)  respectively,  i.e., 


Region  ®  :  r  y 


U  =  -  P(o )f  t>  ('^L'\ntlKp 

*-  27rp  £ro  *W£r0  n«K+2  l-  )  Wk'  ' 


A, 

6(n- - 1 


4  ^e£p(o)f 

zTrr1-  n  to  m2l<*Z  '  '  dA 


‘6  2  ?rr 

Region  ®  :  r  <  R^ 


a  = - * — 

r  2TR,Z 


2.7 rf?/ 


£  p  (o)7  +  Z**P(°)7  -^5-ip(«x-e) 

n To  TTo  '  nsO  h-ZK-ip 

H 

^  /  r  sn-zk-i 

2_  P  (o  /  (~)  p  (<x+Q) 

in'.  ^jL.  n-ZK-|  \f?j  /  'n-Zk  ' 


®  27tP  - 


n-°  K=o  M-°  K*o 


8l« 


2K-I 


(91) 


(95) 


*■  — s  f  *  *p  <«  f  ^  (*r"" 


1  zTrR,2"  «e» 


V" 


2®  Streamf unction 


The  streamf unction  for  the  two  regions  is  found  from  Eqs.  (35)  and 
Eqs*  (9h). 
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Region  ®  :  r  > 


Since  the  differential  equation  in  the  brackets  is  equal  to  zero,  the 
above  equation  becomes 


x A,  K  A1 0 


(102) 


Integrating,  we  get 

(e)  =  \  A,  b0  (co^etC^) 


(103) 


where  the  constant  is  evaluated  from  the  condition  (0)  ■  0 
for  0  ■  0.  Consequently 

<V-i 

And  Eq,  (97)  for  the  st.'eamf unction  in  the  region  r  ;>  R-  becomes 

Jl  ^ 

V,e)s  iA,b0Ke-,)-A,^  |>„4 


(loit) 

(105) 


Region  (b)  :  r  <  R^ 


From  Eqs.  (35)  and  (95)>  we  have 

—  -  —  ZTC^r  UA  Q 
d'r  0 


~  r  a(v  0 

rv. 


ip»i 

0-.0  n  Ki0  £5  n  ^  h-2K- 


iiAJL 


Integrating,  we  get 


t 4^ 


Zp(o)Z  -i-Z**p(o)Z 

ntzk+z  v  ^  h-2K- 


_  ^  A-^er^(o)Z  ^ 


CD 


r_\ 


(106) 


ch*. 


(107) 
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Differentiating  the  above  equation  with  respect  to  0, 

sin^  0  »  1  -  cos*"  0  *  1  -  ,  we  find 

rt  N 


it 

30 


2-  R,1 


/Zvy\Q 


!  i 


oo  x 

tIMPwZ-^ 

htaKtz  ^  r»w^r0  rt-zK-i 


'2.K 


and  substituting 


2CH$ 


-r14LrA-<eZuP(°)£  7— - rf-51-) 

f?  *  nso  n  k -0  (n-2.K+|)(n-2K-iy  \/?  / 


4 


n-iK-l 


(108) 


C^M. 


-C-X)^ 


From  the  second  of  Bqs.  (35)  and  from  the  first  of  Eqs.  (95),  we  have 
that 


at 

a  e 


=  2/)n~2  ar  *1*.  0 


(109) 


-r^A^sI^p.wZ^f— T(^-)  1  *>«£„(<**  8) 


Comparing  Eqs,  (108)  and  (109)  and  rearranging  terms  we  get 


1 £•  *i~4**p-  ‘4  (t-) 


n-iK-i 
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Since  the  differential  equation  in  the  brackets  is  equal  to  0,  the 
above  equation  becomes 

^'(0)  =  o  (1U) 

Integrating,  we  get 

^(®)  =  CS’  (112) 

where  the  constant  is  evaluated  from  the  condition  (0)  -  0  for 
0  ■  0.  Consequently 

3'W  =  C*'°  (113) 

And  Eq.  (107)  for  the  streamfunction  in  the  region  r  <  becomes 


'-K  =  rz  A- ^*9 

‘®  zR,L 


JL 

2 


(11 h) 


Ri  ^  ^  ira  7n  -2Kvl)(n-ZK-i)  V  R,  j  ^ 
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VIII.  SUMMARY  OF  EQUATIONS 


1.  One  Step  in  the  Distribution  of  Sinks 
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SECTION  D 


FLOW  FIELD  FOR  CYLINDRICAL  NON-  UNIFORM  VORTEX  SHEET 

SUMMARY 

In  this  report  the  flow  field  of  a  shrouded  propeller  has  been  inves¬ 
tigated  by  considering  a  vorticity  distribution  on  a  cylindrical  surface. 

For  high  speed  shrouds  the  vorticity  distribution  has  been  assumed  on 

a  cylinder  with  constant  radius  o  which  corresponds  to  the  exit  radius 

o 

of  the  shroud. 

The  vorticity  distribution  on  the  semi -infinite  cylinder  behind  the 
shroud  is  considered  constant  and  on  the  cylindrical  surface  of  the  shroud 
is  approximated  by  a  linear  distribution  and  an  expansion  in  Birnbaum  series. 

The  strength  of  the  constant  vorticity  downstream  from  the  shrouds  is 
calculated  from  the  difference  of  the  exit  velocity  and  forward  velocity  of 
the  shrouded  propeller.  The  strength  of  the  vorticity  distribution  along 
the  shroud  is  calculated  from  the  condition  that  the  shroud  is  a  streamline, 
by  solving  a  system  of  equations  for  the  unknown  coefficients  of  the  terms 
of  the  3irnb--.;m  series. 

For  an  approximate  evaluation  of  the  flow  field  we  might  consider  only 
the  first  three  terms  of  the  series  and  solve  a  system  of  three  equations  for 
the  three  unknown  coefficients. 

Since  the  high  speed  shroud  is  always  a  smooth  curve,  three  flow  conditions 
at  the  leading  edge,  at  the  middle  line  and  at  the  trailing  edge,  might  be 
enough  for  a  good  approximation  of  the  flow  field  by  the  first  three  terms 
of  the  Birnbaum  series. 

The  flow  field,  the  pressure  distribution  and  the  forces  along  the  shroud 
are  determined  from  the  calculated  vorticity  distribution. 
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I.  INTRODUCTION 


In  this  report  the  flow  field  of  a  high  speed  shrouded  propeller  is 

investigated  by  the  method  of  singularity  distributions.  Since  the  high 

speed  shroud  is  a  cylindrical  surface  with  aporoximately  constant  radius,  the 

vorticity  distribution  is  considered  on  a  cylinder  with  radius  p  n  o  = 

o 

constant. 

Making  use  of  non-dimensional  variables,  we  consider  the  shroud  extending 
from  x  »  -  1  to  x  *  1.  Fig.  Dl, 


Figure  Dl 

The  vorticity  distribution  y(x)  on  ‘the  cylinder  o  “  oQ  from  x  =  -  1 
to  x  -  co  is  not  known  and  it  will  be  determined  from  the  boundary  conditions 
and  the  flow  at  x  =  ±  oo  . 

From  the  strength  of  the  vorticity  distribution,  we  determine  the  flow 
field,  the  pressure  distribution  and  the  forces  along  the  shroud. 

The  evaluation  of  the  vortex  distribution  is  facilitated  by  an  approximation 
for  the  representation  of  the  flow  due  to  the  propeller  by  a  semi-infinite 
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vortex  sheet  of  constant  strength  y  B  Y  from  x  =  +  1  to  x  “  +oo  and 

o 

Yq 

a  linear  distribution  y  »  —  (1+  x)  in  the  interval  x  °  -  1  to  X  B  ■+■  1. 


Figure  D2 


In  addition  to  the  above  vortex  sheet,  the  effect  of  the  shroud  is 
represented  by  expanding  the  vorticity  distribution  in  Birnbaum  series 
of  normal  functions  in  the  interval  -1  <.  x  <:  +  1. 

Thus 

**»  =  £  rnN 

n  =M 

where 

*.  <*>  = c, 


YZC*)~  ZzihX1 

Y2  (*)  -  C3  -X  {(-K1 


(2) 


For  an  approximate  evaluation  of  the  flow  field  we  night  consider  only 
the  first  three  terms  of  the  above  series,  and  the  vorticity  distribution 


becomes 


Y(*)=  -^-(n-*)  +  c,fil  t  Cj^Pv'-vCj-jc/TP  I 


U*)~  \  -jx,  \  < 


<  CO 


II.  FLOW  FIELD  OF  A  VORTEX  RING 


(3) 


The  velocity  for  the  flow  of  a  vortex  filament  of  strength  yds  is 
given  by 


->  _  _Y_ 

v  "  ¥7r  RJ 


(W 


where  R  is  the  polar  vector  from  the  vortex  filament  to  the  point  P. 

For  a  vortex  ring  with  radius  pQ  at  the  plane  normal  to  the  x-axis 
at  the  origin,  we  have 


where 


54 


55 


Substituting  Ec.  (8)  into  Eq.  (6),  we  get 
f  aTr 

r  *  ,  *  tr  % 

*4  *  ■'■f  +  (4  -  *■(>  £*  <*+ (fctf) 

Since 


and 


ee»  ‘  ^  5  ^  (%  ~  *?) 

\/  (p>(j  >x)  =  a  ex  +  xr 


we  have  that 

\/*ex=^c  £  V  *  6^,  =  tr 

And  the  velocity  components  in  the  x  and  p  direction  are 


(10) 


(ID 


From  th®  identities  of  Bessel  functions  (Appendix),  we  have 


-six  | 

e  J0  (sP)  J(  (sf.)  5ds  = 


(12) 


-nd 


“  -sj*| 

e  .J(  (sp)  J/s^sds  = 


(n  +  ti  for  K  >  0) 

(»  -  "  for  x  <  0) 


Making  use  of  Eqs.  (12)  and  (13),  the  Eqs.  (10)  and  (11)  for  the  velocity 
components  due  to  a  vortex  ring  of  constant  strength  y  and  radius  p 

o 

at  x  “  0,  become 


.  p  y  r°°  -s|x| 

(('*}-  e  Jo  (*()  J,  (Sf.)  Sots 


/  PJ  f”  -iM 

e  Jt  (sf)  J  (sf,)  sds 


in.  se;i-i?ifi?iite  VDnTF7  SHE^t  4  y  =  y 


For  a  semi-infinite  cylindrical  vortex  sheet  of  constant  radius  o 

o 

and  vorticity  strength  y(x)  ■  y0  in  the  interval  !■£  x<e©,  the  flow 
field  is  found  by  integrating  Ecs.  (11;)  and  (l£)  for  u  and  v  from  1 
to  <=o  with  resoect  to  X. 


Thus, 


-5|X-X'| 

e  (Sf)  J  (sf.)srf  Sol*' 


f (Xj  p  y  600  —  sjy  — ,x/| 

vCf’*)  ‘J  t  j  e  J,  (s(f)J,  (Sfo)  Sets 


("+"  Tor  (x-x1)  >0,  for  (x-x')<0) 

From  the  evaluation  of  the  integral  with  respect  to  x'  we  have  (Ref .Dl) . 

I  -s(X-|)  /? 

f,CC  -.s|X-X'|  y  5  £  ^ 


P  -«s|X-X'|  ;  /  S  C 

I  V  -  — /e'S^"-[V±i])  ^ ' 

s  ^  U-y »_» 

Substituting  into  Eqs.  (16)  and  (17),  we  find 


'Z  <.* 
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/co 

e  J,  (Sf)  T,  (spoys 

'o 

M.  (p,x)=  -  ^j[  [e  SC  -2  T0  J,  (s£o)cfc  jb  X  *  | 

id 

V  ((Yx)=  -  j  e  ^  J,  (s(f)  J,  <Sfo)  c^s  ^  ^  I 

/  A 


(q*): 


.&if’.-s(’,-,> 


J|  (£f)  Jj  tefo)  «ts 


From  the  above  Eqs.  (19)  and  (20),  we  see  that  for  x-*tco  the  velocity 
becomes 

r  co 

«.«.  =  e  J0  (60  (sft  =  0 

'o 

/ 

-CO 

"tco  =  '  Z  e  l  J'l  ds  +  f»  ^  J.  (Sf)  3J  (5?o) 

'  J,  L 


+  fo  X> 


e-c< 


f- 


IV.  FINITE  TRLWTLAF  VOITEX  SHEET  *  -  h-  Q  +  * 


For  a  finite  vortex  sheet  of  constant  radius  pQ  and  strength  tf(x)= 
in  the  interval  -1  <.  x  <  1,  the  flow  field  is  found  by  integrating  Eqs,  (Hi) 


and  (15)  for  u  and  v  from  -1  to  fl  with  respect  to  x. 
Thus, 

f+lpy  ( 10  “Slx-'x'l 

«(f I *)  -/  V  0+  X')  e  J0  (Sf)  J  (sf„)  sd  so/  *' 


/  ,  f*1  P)  f®  -s|x-*'| 

±  I  (i+*')e  J(  (sp)j;(s{Osols«l  x' 

I  4  ("+  "  for  x  -  x1  >  0,  for  x  -  x1  <  0) 


From  the  evaluation  of  the  integrals  with  reject  to  x1  we  have  (Ref.  Dl) 


2 

O 


x  ^ 


-s|x-x'| 


d  C04- > 


j£S^i.SX  ->(“) 

>  -1^  x<  tl  (25) 


Ji  -SxJ  j  (v) 


s 


‘  I  ^ 


|z  “  [j  sj |  jti. 


X  -| 


'+i  -sjx-ix'/  rr[*-e  + 

x'e  cty'=  <  (26) 

I  •f't’T ‘t'e'*5(,tT)c<^"^sXJ  -+(v) 


2.  [ockl.^5  -  -~-xu«^ls 
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Substituting  into  Eqs.  (23)  and  (2h)  for  u  and  v  we  find 

Ado  p 

e  cc&Js  J0(sfi  J(  (s^0)M,  j 


eta-Ti  ■  *  +  '-« 


j(l  t- 3-)  sx  -i-  axj  sx  j  J0(sf) ^  (sp)ct 


P  Vo  -STCf  1 

_l  e  (l-  -~)a^as  +•  C&lJIs  J6(^f)  J,  (Sf>0)^S 
/o  L  J 


eA  ** 


^  =  e 


-if  . 


(n-~)A^is'C^Js  (Sf)  Jj  (S^0)  c=4s  %<-| 


00  r 

-Sr  n  3 

■6  (l  AX  ~  i  3j  Jj  ^ 

;0  L  J 

f' 


'kx<H 


_5-Xf  -1 

£  (I-  /U/xx.  I  S  tCX^Js  J(  ($£>)  J(  (5^  ^  /< 


V.  FINITE  VORTEX  SHEET  <  Ux)  =  C 


For  a  finite  vortex  sheet  of  constant  radius  p  and  strength 

0 


in  the  interval  -!•£  x^  +  |  the  flow  field  is  found  by 


integrating  Eqs.  (Ill)  and  (l£)  for  u  and  v  from  -1  to  +1  with 
resuect  to  x. 


Thus, 


Tt  (sf)  J,  (sfc)  s  cfc  oCf! 
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+7  f  W  . 

SV  I  7T5T- 


/  ±  2. 


I-*' 

I 


&  J,  (Sf)  X,  (Vo)  s  di  dn' 


(30) 


("+  n  for  x  -  x*  ^>0,  "-M  for  (x  -  x1)  <  0) 

From  the  evaluation  of  the  integral  with  respect  to  x‘  we  have,  (Ref.Dl) 


■ti 


^5T-&-slx-*W  = 

/  +  *'  e  ^ 


-l 


Or  r-  /  *\ 

Ta  Fa(-S) 


A 


y< -/ 


-s* 


jl,  Fn(s)0+fT+F„«(! 


(31) 


±  a  *[$  I  V)0tf)'"  t  F.  V)  3^"  f] 


I 


ft* 


F0(A) 


wnere 


/+* 


fv^w] 

(v) 


-jifi'  I  <  X 


p  .  rfn  a(h+fe)* 
1  I  ^TTTa 


F"  ^j  =  ifc 


i+"2 


L  -v*c/ 


A-i2 


(32) 


F(°0  =  4£aJl  -iV 
n  2^  +  / 


«,<•>-  if 


Z  o* 


Substitutin';  irto  Fqs.  (29)  and  (30)  for  v  and  v  the  velocity  field 
becomes 


r00 


a  =  A-L-.I  (-5)  J0Csf)T  (Sfo)  Sds 


X  <  - 


6l 


k«£l 

2. 


+  e  \  Z  Fn  (-s)(i<-t1)  n]|- 


■J0(sf)  J,  (s&  eds 


e.s  -s* 


£  *"  FI  (s)  (sf)  J  (s$)  ids 


Is  -x 


e*  V  F"  (-S)  T  (Sf)  J,  (sf.)  icU 


f>  C. 

v=i^i 


“ £SZ[fo("'s)&^',S  +■  t  t (-s)0+fz)  ]]• 


I  («f)T  (Sfo)  Sdi 


-IfXSl 


£  7r  Fo  (s)  J(  6^)  jo  Csf0)  sds 


/  5  X 


rPTTK  TIPTKX  PHgm7y(x)  =  c  ^\h-vl 


For  a  finite  vortex  sheet  of  constant  radius  oq  and  strength 
v(x)«  CxyT-?in  the  interval  ,  the  flow  field  is  found  by 

integrating  Eqs.  (Hi)  and  (15)  for  u  and  v  .  from  -1  to  +■  1  with 


>C  f°°  I - 72.  -Sl'K-X'l 

IT  J,  (*()!,(*{.)  !<*•**? 

/  0 


«r.*r, 


R5 

("+  "  for  x  -  x'  ?  0,  for  x  -  xT  <  0) 


From  the  evaluation  of  the  integral  with  respect  to  x>  we  have  (Ref.  Dl) 


TrA  (-s)  e' 


x  <  -j 


I - 7  -S  X-X'j  , 

/  ^/-x/2  a  <3^  = 


|e  -h Bnte)*h 


i'll  ,  it  f\ 


+  "  for  x  -  x»  >0,  for  x  -  x'<  0) 


7r  A  (s) 


jU,  -/<%<  i 


I  <,  x 


where 


L 


I  g2-1^  Brnm 

2-  A7|  IST2.  oo  +  l 


05  o-i  1-2  K 

(n-11-iK)!’  VTTkTT 

2.  K 


00  o0  +  Zfc+j  | 

B  (a)  =■  2  t — -.j-- -  A  fn  z  even) 

n  K*o  (n+lk+l)!  n+ZK*<-3  -j+K/  ^  ' 
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1 S  '  I  -  f-  x.A 


r»)*  ^  U^-^-hj 


V£  ^  fK’“  V 

Substituting  into  Eqs.  (35)  and  (36)  for  u  and  v  the  velocity  field  becomes 


pc. 


«s  x 

e  r  A  (-s)  J„  (sf)  T,  (s^)  sds 


x<  -I 


-s*r  _  ,  N  v  3  1 

■  e  A(s)(£+^  £  6„ (s) Tt11 

3  L  L.  J 

+  eSX[A(-s)(? -^'•x -a  f^)H> -*ikZ  B„ 6s)*n]J . 


kI0(Sf)  Kspjsds 


"S* 


e  'Jr  A Cs)  J0  Csf)  J  (sf0)sc& 


TT  a  -  I  e 


** 


^  A  (-s)  I  Csf)  Jo  (sf0)  sdi 


X<  - 


V=  ~J  l^jACs)^^'1  *  +  xf^)-(l-H*f/lZ  6rt(s)xn] 

-  e  A(-s)(~l-/uv,  ’y  -  % ^/- -x X) f  (/ -  OC2-)3^  Z.  8h(-s)*hjj* 


^  (sf)  (Sfo)  Sols 


-!<•%< 


(ho) 

!<>%<  + 1 


"’^'X 

e  9r  Acs)  Jjfsp  Jj  (s£>0)  scU 


/<  * 
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where 


a. 


6xk 

Ktl 


tad  Ba,k  Vk  are  clven  b^r  Ef,s-  (38). 

Substituting  into  Bps.  (1,1)  and  (1,2)  for  u  and  v  the  velocity  field 
becomes 


U 


P(s)  Jg  Of)  Jj  (sf\>)  sds 
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And 


rP(-s)  J,  (Sf)  J,  (sfo)sds 


TC  -I 


e^TpwC-r  -tv^'h i  QnCs)'Vr‘] 

-  e  *[p<-s)(f-  ^V*Vmc*)+  (i-^2  C-s)^h]j*  (w) 

•J,  (sf)  j;  (5fo)  sots 

e'5*?'-  PCs) Jj (sf)  Jj  (s(>o)sols  -^*61 

yjTT  M^pJ_Tnn"r?APy  ypnnvRy  SH^TT 

Ter  the  representation  of  the  flow  field  cf  a  shrouded  propeller  by 
singularity  distributions  we  determine  the  strength  of  the  vortices  fro:.;  th° 
fact  tint  the  shroud  should  be  a  streamline  of  the  flow. 

Si r.c'-  only  high  sp'-'ed  shrouds  are  considered  in  this  paper,  the  vorticit; 
distribution  has  been  assured  on  a  cylindrical  surface  of  constant  radius 
which  corresponds  to  the  exit  cylindrical  radius  of  the  flow. 

The  strength  of  the  semi -infinite  vortex  sheet  v  is  evaluated  fro n 
the  forward  and  exit  velocity  u  and  u  ,  i.e., 

(J  fc* 

\  -  “e  -  =  u+ „  -  «  „  «  *>e  " 

The  coefficients  Cp  c^,  and  are  evaluated  from  the  satisfaction  of  the 
boundary  conditions  on  the  shroud.  Thus,  for  a  number  n  of  coefficients,  we 
need  a  system  of  n  boundary  conditions. 
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I 


In  our  special  case  with  n  =  3  we  need  three  boundary  conditions. 


It  might  be  recommended  to  take  the  conditions  at^s-^  (-i)J 
D**0 )  f  (°)]  )  and  [x »i,  ^  3^»s  (i)J  with  the  corresponding  slopes. 

oif  (**&)  _  v  (*,  el 

ot}L  ~  u  (*) 

where  the  velocities  v(x,p)  and  u(x,p)  are  given  by 

=  u,  (%,£)  t  6Cj,C?(,f)  +  n3 (*,(>) *■“(*>(>)+ 


v  (*>f)  =  v,  (x,f  )  f  ^  (*,f)  t^3  (**f)+/i/  (*>f )  1-  W*,f) 

•  0(  f 


Thus,  we  have 

a 

} 

cot 

%  -~l  1 

f  =  fs  <*0 

V 

X  =  -°<o 

cut 

X  =  0  $ 

'  £  *  f  4  («) 

:£.=-*  -o 

U  +-|  ~ 

cot 

X-sf  |  ^ 

'  C  *  &  (♦') 

Hairing  use  of  Eo s.  (1*8)  and  (i;?),  Eo.  (ii7)  becomes 


(U7) 


(18) 


(1 h9 ) 


(^CK  +  v$)  t  C,  V,  f  cz  +  c3  V3  1 
K  +  +■  uq)  +  C.|  u,  TCZUZ+  Cj  U3 


CK 


**°>  f  sf.w 


r  ^g)  f-  0/  Vj_  +-  2^2  *  ^3  ^3 

( wo^  ^  +*  wp)  +  Cj  U|  u-  cz  c<  j,  +  C3  C-t3 


L 


•  ,-«,  =  0 


(50) 


Solving  the  system  of  Eqs.  (50),  we  find  the  values  of  c  ,  c  and 

1  2 

which  determine  approximately  the  flow  field  of  the  high  speed  shrouded 
propeller. 


c 


3 
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IX.  APPENDIX 


IDENTITIES  OF  BESSEL  FUNCTIONS 


From  the  identity  of  the  Bessel  functions  we  have 


-  I  e  S  Jo  ds 


a  >  o 


Differentiating  with  respect  to  a,  we  find 


d  i  _  __  a 

s «  fiFiT*  ~~  C^Vk1)' 


00  -so. 

e  T  (sb)  sots 


Substituting  a  by  x  and  b  by  r,  we  have 


L-1**1]* 


/  -sm 

=  ±  I  e  1  (sr)  sds 


where  the  "  si,  r.  is  for  x  p»  0  and  the  sign  for  x 0 


r7  =  ^  d 


Multiplying  Eg.  (53)  by  cos  91  end  integrating  from  0  to  2tt  with  respect 
to  0 1  we  ,  ■  t 

r27r  ni  ftrrfc D  .S|X1 

j  ^9'=  t j  j  e  Jo(sr)sds[aw6']^6'  (55) 

Making  use  of  the  following  identity  for  the  Bessel  functions 

J  <»  =  J0  (sf)  l  ( sfo)  t  2.  £  Jn  (sf)  "  6'  (56) 


and  of  the  orthogonality  property  of  the  trigonometric  functions,  we  perform 
the  integration  with  respect  to  O’  and  we  find 

fz^~  i  r°°  —  s  i*xi 

?  ww**  (S7) 
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Similarly,  substituting  a  .by  .  |x|  and  b  by  r  in  Eq.  (£l),  we  have 

^  “'’I*1  r  /  v  i 

e  I0  C^Oots  (58) 


Differentiating  Eq.  (£8)  :4th  respect  to  p  ,  we  find 


3e°  i/TsF  ~"'1  ^hTW-"] 


(5?) 


Integrating  from  0  to  2tt  with  respect  to  0*  we  get 


„  A,  ,  fair  r oo  .  . 

e0-e^&  ;V. 

FF  1  .  -f 


['•H 


t?6s 


(60) 


Substituting  Eq.  (£6)  and  interchanging  the  integration  and  differentiation 
in  the  right  member  of  the  above  equation,  we  perform  the  integration  with 
respect  to  91  and  making  us?  of  the  orthogonality  property  of  the  trigonometric 
functions,  we  find 


Performing  the  differentiation  with  respect  to  p  ,  we  get 

o 


(61) 


ol  r 


f- 


J0  (sf»)  =  -  si,  (Sfo) 


(62) 


Substituting  Eo.  (62)  into  Eo.  (61),  we  find 


1  •  /  e0-e  ^  6 


s/*|, 


2.7r 


(63) 
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SECTION  E 


VELOCITY  FIELD  AND  FORCES  ON  A  FLAT  PLATE 
DUE  TO  A  CYLINDRICAL  JET 

I.  Elliptical  Cylinder  With  Major  Axis  Parallel  to  the  Flow 

In  this  section  the  pressure  is  found  on  a  plate  parallel  to  the 
flow  and  normal  to  a  cylindrical  jet  with  elliptical  cross  section. 

Since  the  velocity  of  the  jet  is  very  high  in  comparison  to  that  of 
the  flow  field  around  it,  we  consider  that  the  jet  is  a  rigid  cylindrical  body. 

The  potential  flow  field  around  a  cylinder  with  an  elliptical  cross 
section  is  calculated  from  the  theory  of  conformal  mapping,  making  use  of 
the  complex  potential  for  the  flow  field. 

For  an  ellipse  in  the  z  (x,  y)  plane  with  major  semi-axis  a  and 
minor  semi-axis  b  the  transformation 


maps  the  ellipse  into  a  circle  with  radius  R  ■  j  (a b ) 


*The  complex  potential  for  the  flow  velocity  in  the  'j  -plan®  is  given  by 


F(0’v(r+-f) 

where  ^  is  given  by  Equation  (1). 

Solving  Equation  (1)  for  J  we  find 


r- 


5  +  -  c2- 


where 


"L  ^  i2- 

rr  Dl  —  b 


(4) 

(5) 


Substituting  from  Equations  (4)  and  (5)  into  Equation  (3)  we  find 


Making  use  of  the  elliptical  transformation 


2  =  -  c  a (£+L7)  <7> 

we  get 

#  =  C  <x>*4  J  <**7  ;  (8) 


Substituting  Equation  (7)  into  Equation  (6)  the  complex  potential  for  the 
flow  velocity  becomes 


Performing  the  algebraic  manipulations  and  making  use  of  the  identities 
of  the  hyperbolic  functions,  the  above  Equation  (9)  becomes 


Separating  the  real  and  imaginary  part  of  the  complex  potential  we  find 


Thus,  the  velocity  potential  and  the  stream  function  in  elliptical  coordinates 
are  given  respectively  by 


The  velocity  components  become 
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where 


(16) 


Since  the  flow  is  isoenergetic,  the  Bernoulli's  equation  gives 


c  o'n  s 


-t 


(17) 


where 


7. 


(18) 


The  pressure  field  on  the  plane  normal  to  the  elliptical  cylinder  is  given  by 


(19) 


Making  use  of  Equations  (14),  (15),  and  (18)  the  above  Equation  (19)  becomes 


The  force  on  the  plate  is  found  by  integrating  the  pressure  field  given 
by  the  above  equation  all  over  the  plate  outside  the  ellipse. 


For  an  elliptical  plate  with  the  semi-foci  distance  c  as  that  of  the 
jet  and  with  major  semi-axis  A  and  minor  semi-axis  B  the  force  is  given  by 
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Substituting  for  &&  from  Equation  (20)  we  get 


Expanding  the  integrand  in  Equation  (23)  and  making  use  of  the  identity 


Equation  (23)  becomes: 


Integrating  with  respect  to  we  find 
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Since  we  have  that: 

*  ~  STfc  > 

-\{Ot+b/ 

4 


the  above  Equation  (26)  yields 


L  =  -{v 


JC 


~2ti  4-B 


A-t-B 


Oc 


-i> 


OL-l-b 


— 


(27) 


(28) 


And  the  lift  coefficient  becomes 


where  S  is  the  reference  wing  area 

$  =  TtfAZ-otij 


* 

Taking  as  reference  area  S  the  area  of  the  elliptical  cylinder, 

* 


where  S  is  given  by 

"  JVot  b 
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II.  Elliptical  Cylinder  With  Major  Axis  Normal  to  the  Flow 


For  an  ellipse  in  the  z  (x,  y)  plate  with  major  semi-axis  a  and 


minor  semi-axis  b  ,  the  transformation  given  by  Equation  (1)  maps  the 

1 

ellipse  into  a  circle  with  radius  R  B  %  (a-/-b) 


A  V  fu- 


Making  one  more  transformation  given  by  Equation  (31) 


we  find  that  the  velocity  U  is  turned  by  90° 


where 


are  given  by  Equations  (1)  and  (31). 
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Substituting  from  Equations  (4),  into  Equation  (32)  we  find 


Making  use  of  the  elliptical  transformation  given  by  Equation  (7) 
and  performing  the  algebraic  manipulations,  the  above  Equation  (33) 
becomes 


Separating  the  real  and  imaginary  part  of  the  complex  potential  we  find 
the  velocity  potential  and  the  stream  function  in  elliptical  coordinates 
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Making  use  of  Equations  (37),  (38),  and  (18),  Equation  (19)  becomes 


Performing  the  integrations  we  find  Equation  (26). 

Thus,  the  force  on  the  elliptical  plate  and  the  lift  coefficient  are  given 
by  the  same  Equations  (28)  and  (29)  as  in  Section  EL 

For  major  semi-axis  b  and  minor  semi-axis  a  Equations  (28), 
(29),  and  (30)  become  respectively: 
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SECTION  F 


EXPERIMENTAL  DATA  ON 
JET  INTERFERENCE  EFFECTS  ON  VTOL  AIRCRAFT 

The  jet  issuing  from  a  VTOL  aircraft  such  as  the  fan- in-wing 
type,  induces  pressures  on  the  lower  surfaces  of  the  aircraft  which  in 
most  cases  cause  a  loss  in  lift  and  a  nose-up  pitching  moment  during 
transition.  Simple  experiments  designed  to  study  the  effect  of  a  jet 
emerging  normal  to  a  flat  plate  placed  parallel  to  an  air  stream,  and 
similar  tests  on  a  solid  circular  cylinder  normal  to  the  plate  were 
carried  out  at  the  71  x  10'  tunnel  at  Langley  Field.  These  data  have  not 
been  published  and  were  obtained  through  a  visit  by  Republic  personnel. 
The  pressure  distributions  found  on  the  plate  are  shown  in  FiguresFl  and 
12  for  the  jet  and  the  cylinders,  respectively. 

Positive  pressures  are  produced  ahead  of  the  cylinder  with  rather 
high  negative  pressures  to  the  side  and  rear.  Similar  results  may  be 
seen  for  the  jet  with  the  exception  that  the  negative  pressures  are  higher 
and  extend  to  a  greater  distance  from  the  center  of  the  jet.  The  result  is 
a  predominance  of  negative  pressure  to  the  rear  of  the  jet  which  in  the 

case  of  a  VTOL  aircraft  in  transition  would  produce  a  loss  in  lift  and  an 

/ 

increase  in  nose-up  pitching  moment.  Quantitative  results  for  the  two 
pressure  distributions  are  given  in  TableH  in  coefficient  form  iof  the 
front  and  rear  halves  of  the  plate  area. 
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In  the  case  of  the  cylinder  the  lift  on  the  front  half  of  the  plate  (Lj) 
is  negative  (C^  -  anc^  acts  at  a  P°^nt  radii  downstream  from 

the  cylinder  center.  This  checks  reasonably  with  potential  flow  theory 
which  indicates  a  lift  coefficient  of  somewhat  below  -1.  The  lift  on  the 
rear  half  of  the  plate  (L^)  is  much  more  negative  (Cj^  a  -5.59)  and  the 
center  of  pressure  of  this  lift  is  at  a  point  about  two  radii  behind  the  center. 
The  total  result  is  a  negative  lift  of  =  -3.  1  acting  at  a  point  two  radii 
behind  the  center,  thus  producing  a  nose-up  pitching  moment. 

The  case  of  the  jet  is  similar  with  the  exception  that  higher  values 
of  the  negative  pressure  peaks  at  the  sides  of  the  jet  predominate  to  the 
extent  that  the  pitching  moment  of  the  front  half  of  the  area  (Mj)  in  this 
case  becomes  negative.  Also,  the  total  negative  lift  is  about  twice  that  of 
the  cylinder  and  the  total  force  acts  at  a  point  about  12/3  radii  behind  the 
center,  producing  also  a  larger  nose-up  pitching  moment. 

The  higher  negative  pressures  induced  by  the  jet  as  compared  to 
the  case  of  the  solid  cylinder  is  attributed  to  a  downward  "pumping"  action 
in  the  mixing  region  of  the  jet.  Measurements  taken  at  static  condition 
indicate  a  total  thrust  loss  on  the  plate  of  only  one  percent  of  the  jet  thrust. 
The  pumping  action  of  the  jet  increases  with  forward  speed  due  to  the  mixing 
action  which  is  also  responsible  for  the  drag  on  the  jet.  The  coefficients 
CL  and  are  obviously  related  to  the  drag  coefficient  of  the  jet 

as  they  are  all  produced  by  a  negative  region  behind  the  jet.  The  effect  of 
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Reynolds  number  may  be  inferred  from  Figure  F3  where  the  pressure 

distribution  on  the  plate  near  the  cylinder  is  compared  with  that  on  a 

cylinder  taken  from  Reference  FI.  The  drag  coefficients  (Cj-j«il.3 

based  on  cylinder  length  and  diameter)  obtained  are  characteristic  of 

a  laminar  flow  separation.  This  is  to  be  expected  since  the  Reynolds 

5 

number  is  below  the  critical  value  of  about  3.  6  x  10  (based  on  forward 
velocity  and  cylinder  diameter).  For  most  VTOL  applications,  the 
Reynolds  number  will  be  above  this  value  and  consequently  a  turbulent 
separation  will  produce  a  smaller  area  of  negative  pressure  behind  the 
cylinder.  If  this  reasoning  is  applied  to  the  case  of  the  jet,  smaller  lift 
losses  and  reduced  nose-up  pitching  moments  are  to  be  expected  with 
larger  forward  speeds. 

In  conclusion,  the  values  of  lift  and  moment  coefficients  determined 
from  integration  of  the  pressure  distributions  shown  in  this  report  and 
the  lift  and  moment  coefficients  from  the  momentum  theory  outlined  in 
Reference  F2  have  been  computed  and  the  results  are  presented  in  Figures 
F4  and  F5.  Also,  some  experimental  values  from  Reference  F3  are  shown. 
The  losses  in  lift  and  the  increase  in  nose-up  pitching  moment  caused  by 
the  jet  interference  are  sizeable  and  could  easily  explain  the  actual  effects 
observed. 

The  conclusions  concerning  jet  interference  which  may  be  drawn 
from  this  study  are: 
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The  jet  issuing  from  a  VTOL  may  be  regarded  as  a 
solid  cylinder  with  a  downward  pumping  action  in  the 
mixing  region. 

The  flow  on  the  front  half  of  a  plate  normal  to  a 
cylinder  resembles  a  two-dimensional  potential  flow 
about  a  cylinder. 

The  losses  in  lift  and  increase  in  nose-up  pitching 
moment  caused  by  the  interference  of  a  jet  are  large 
and  do  explain  the  deviations  from  momentum  theory 
observed  in  experiments. 
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SECTION  G 


DESIGN  OF  DEFLECTION  VANES 

The  design  of  a  grid  of  turning  vanes  constitutes  a  part  of 
the  general  VTO  problem.  In  particular,  it  is  understood  that  it 
is  entirely  undesirable  to  expose  a  propeller  or  fan  to  a  flow  at  an 
angle  with  the  axis  even  for  a  short  period.  A  system  of  baffles 
or  turning  vanes  may,  therefore,  be  required  in  any  practical 
solution  of  the  VTO  problem.  We  shall  not,  in  the  present  paper, 
specify  when  or  where  such  a  baffle  system  is  or  may  be  required 
but  merely  give  the  design  information  required  for  the  proper 
solution.  The  proper  use  of  deflection  vanes  is  taken  up  in  Section  I. 

The  basic  conformal  transformation  for  cascades  has  been 
known  for  some  time.  (See  Ref,  Gl.  )  We  shall  in  the  following  give 
the  direct  routine  method  for  the  solution  of  any  given  case  based  on 
methods  given  in  Ref.  Gl. 
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In  Figure  G1  is  shown  a  set  of  turning  vanes.  Inlet  velocity 
is  Wj  at  an  angle  shown  as  OC  ^  and  an  outlet  velocity  w^  at  an 
angle  .  in  general,  it  is  desirable  to  provide  at  least  a  small 
acceleration  of  the  flow  through  the  vanes,  which  means  that  w  is 

Ld 

to  be  a  little  greater  than  w^  or  the  angle  0(  is  to  be  a  little 
larger  than  OC ^  •  This  also  means  that  there  is  a  small  overpressure 
on  the  upstream  side  of  the  grid.  In  fact 

f, -fx  -  jf  ({ wi~  w?) 

This  may  also  be  written 

where  w^  is  the  tangential  component  of  w^  and  w^  [ s  the 
tangential  component  of  w^  .  (See  Figure  G2.  )  The  resulting  pressure 
force  as  can  be  seen  is  normal  to  the  grid.  The  force  along  the  grid 
axis  is  further 

where  wn  is  the  velocity  perpendicular  to  the  grid  axis  and  common 
to  inlet  and  outlet,, 

The  right-hand  section  of  expression  (1)  can  be  written: 

)(W*t  +  W'*) 
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The  geometric  addition  of  the  forces  (1)  and  (Z)  therefore  is: 

^>(^21:  -f  Wa^.  -f  Wit 

From  Figure  G2  it  may,  however,  be  seen  that  the  bracket  expresses 
the  mean  value  of  the  velocities  Wj  and  w^  .  Let  the  vector 

■=■  wj  f  Wa. 

Then  the  force  on  the  grid  for  a  single  element  of  width  t  and  unit 
length  is 

-  A  .  .  ,  . 

(3) 


F=  (wlt-wn)vjw 


This  force  is  seen  to  be  perpendicular  to  the  mean  or  average  velocity 
wm  as  shown  in  Figure  G2„ 

It  is  interesting  to  note  that  the  circulation  ["*  can  be  read  off 
Figure  G1  as 

r-  ^ )t  (4) 

This  is  true  since  the  elements  cancel  everywhere  except  at  plus  or 
minus  infinity,,  Note  that  the  force  p  fn  Figure  G2  points  downward. 
Hence  a  counter  clockwise  circulation  P 

Thus,  finally  P  =•  yO  | 


P  ==•  yO  P  W yy,  and  P  -L 


(5) 
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with  Vv/y^  —  ~h 

It  is  at  this  point  interesting  to  make  the  remark  that  this 
result  is  true  regardless  of  the  shape  of  the  grid  elements.  Since 
the  force  is  perpendicular  to  the  mean  velocity  the  center  of  gravity 
of  the  vortex  system  is  located  behind  the  semichord  of  the  elements 
in  FigureGh  However,  the  exact  location  is  not  of  importance  as  is 
the  case  for  a  single  wing  or  wing  section. 


It  is  desirable  to  know  the  value  of  the  lift  coefficient  C 
since  experience  shows  that  excessive  values  are  not  allowable, 
may  write: 


L 

We 


(6) 


where  X  is  the  chord  of  the  element  (Figured).  In  contrast  to  the 
expression  (5)  which  is  an  exact  expression  regardless  of  the  shape 
of  the  element,  the  value  of  as  defined  by  Equation  (6)  is  strictly 
accurate  only  for  small  deflections  and  infinitely  thin  elements. 


In  Figured  it  is  seen  that  the  value  w^  is  very  much  smaller 

than  the  velocities  w.  and  w_  .  In  actual  practise  the  value  of  w 

1  2  r  m 

would  be  made  slightly  greater  than  Wj  and  smaller  than  w^  to 

provide  continuous  acceleration  through  the  channel.  The  expression 

2 

(6)  where  is  based  on  w^  would*  therefore,  give  the  wrong 


impression  of  the  realistic  value  of  C  .  Since  nevertheless  we  shall, 

L 
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for  convenience,  adopt  the  use  of  Equation  (6)  we  shall  remember 
that  the  real  value  of  C  is  very  much  lower  for  large  angles  of 

Lj 

deflection.  In  fact,  the  velocity  in  the  representative  point  may  be 
defined  as 


W,  j  +  |  W2 

5 


and  one  may  write  with  the  realistic  value  of  C 


p  = 

4, 

w 

!  and  | 

w„ 

1  1 

2 

W,  +  J  Vv/2/ 


(?) 


represents  the  numerical  values  of  the 


vectors  w^  and 


Consequently,  the  "realistic"  value  of  the  lift  coefficient  is 
in  terms  of  the  fictional  value 


w,  +  \a/2 

jw,J  + 1  w2j 


(8) 


The  ratio  of  the  real  to  the  fictional  value  is  thus  given  by  the 
square  of  ratio  of  the  geometric  sum  of  the  vectors  w^  and  w^  to 
the  square  of  the  arithmetic  sum  of  same:  Since  the  angles  (X  and 
(X  in  the  cases  of  importance  in  VTO  design  are  of  approximately 
the  same  value,  we  shall  employ  the  mean  or  average  value  OC 

m 

We  have  then  approximately 
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where  thus  (  2L—  °(rr\ )  is  one  half  °f  the  total  deflection  angle.  In 
Z~ 

o 

the  two  cases  shown  in  the  following  the  deflection  angles  are  27 

and  55°  and  the  "real"  values  of  C  are  therefore  lower  than  the 

Xj 

2  o  2  o 

given  artificial  values  by  factors  cos  13.5  and  cos  27.5  or 
0.94  and  0.79,  respectively.  In  the  second  case  shown  there  is  an 
indicated  value  of  C  -  1.  04,  the  "realistic"  value  with  uniformily 

Xj 

contracting  channels  is  therefore  only 


C  =  1.  04  *  0.  79  =  0.  82 


Reverting  to  the  artificially  defined  value  of  C  (Equation  (6))  and 

L 

using  Equation  (5)  one  has  the  relation  for  C  (without  correction) 

■Lj 


c  Z  Z  t  ~  Wfc 


l 


lA A 


Yn 


EXAMPLES 


As  experimental  checking  is  proposed.we  shall  show  two 
examples  obtained  by  routine  method  from  Reference  Gl.  Figure  G3 
shows  a  grid  for  a  total  deflection  of  27°.  The  angles  with  the  normal 
are  chosen  to  be  10  on  the  inlet  side  and  17  on  the  outlet  side,,  The 
spacing  t  is  taken  as  30  m/m  and  the  projection  of  the  chord 
^  -  40  m/m.  We  are  performing  the  routine  calculation  for  circular 
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arc  profiles.  The  values  of  w  ,  w^ ,  and  w^t  may  now  be 
immediately  taken  from  a  graph  similar  to  that  of  FigureG2and  we 
have  for  C 

L 

CT.  a  0.72 

l  uC 

The  only  other  value  required  is  the  radius  of  curvature  of 
the  circular  arc  section.  This  is  obtained  by  using  Figure G4 which 
is  based  directly  on  a  routine  method  given  in  Reference  Cl. 

One  has  only  need  for  two  quantities.  One  is  the  mean 
deflection  angle 

/,=  3.5° 


The  second  is  a  factor  which  contains  all  the  complications  of 
the  theory  but  which  may  be  given  in  a  single  chart,  Figure  G5,  taken 
from  FjferenceGl.  With  s  3.5  and  jL  a  0.75  one  gets 

directly  from  this  graph  yU.  -  0.675.  This  value  is  now  used  to 
find  the  angle  "exaggeration"  up  stream  and  down  stream  by  the 
following  scheme  (see  Figure  G4) 
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A»-/e 


3.  Sa-HO°)  =  2C>- 


Thus 


o 

23.  5 


The  mean  radius  of  curvature  is,  therefore,  (see  Figure  G3) 


R= 

7  s/n  2o° 


58.  5  m/ m 


This  is  the  value  for  the  radius  of  curvature  of  a  skeleton  circular  arc 
element.  In  reality  the  channel  will  be  designed  for  a  gradual  increase 
in  the  velocity,  thus  being  slightly  thicker  in  the  middle  section.  The 
total  deflection  designed  for  27°  will  be  somewhat  smaller  as  a  result 
of  the  effect  of  a  boundary  layer  on  the  upper  rear  surface.  The 
projected  value  of  C  =  0.72  (to  be  corrected  by  a  factor  of  0.94 
as  above)  will  be  somewhat  decreased  in  an  actual  experiment. 

Finally,  a  second  example  has  been  worked  through  also  based 
on  the  ro'  tine  method  of  Reference  Gl.  The  total  deflection  is  here  55° 
(see  Figure  G6).  The  projected  chord  length  is  again  40  m/m  and  the 
spacing  is  chosen  as  20  m/m  ,  Repeating  the  same  calculations  the 
results  obtained  are  shown  in  this  figure.  The  resulting  radius  of 
curvature  is  33.6  m/m.  Experiments  will  be  proposed  to  confirm 
the  accuracy  of  the  calculations  and  to  obtain  knowledge  of  the  optimum 
efficiency  possible  by  refinement  in  the  design  procedure. 
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Experimental 

Investigation  of  a  cascade  of  circular  arc  sections  with  and 

without  a  downstream  channel 

[Project  No.  6-C  3338  -  548722  of  Republic  Aviation  Corporation 
Fariningdale ,  New  York,  U.S.A.,  Jan.  18,  1 96 Ijf 
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Summary : 

A  cascade  of  cirgular  arc  sections,  which  shall  serve  as  a  deflect¬ 
ion  grid  in  a  55  -bend,  was  investigated  by  the  wake  traversing . 
method.  Measurements  were  taken  with  and  without  a  channel  down¬ 
stream  of  the  cascade. 

It  was  found  that,  though  the  measured  discharge  angle  fitted  very 
well  with  theory,  the  losses  were  rather  high.  They  still  increased 
when  a  channel  was  installed  behind  the  cascade.  The  reason  is  that 
separation  occurred  at  the  lower  surfaces  (pressure  sides)  of  the 
blades.  V/hen  using  a  downstream  channel,  the  inlet  air  angle  of  the 
cascade  decreased  and  separation  increased. 
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1 )  Introduction 


Deflection  vanes  are  commonly  used  in  bent  channels  to  avoid 

* 

wall  separation  and  thus  to  reduce  total  head  losses.  For 
reasons  of  convenience,  often  blades  of  circular  arc  sections  : 
arc  used.  This  report  presents  measurements  with  a  cascade  con-  j 
sisting  of  such  blades  which  is  to  be  used  in  a  35°-bend.  The 
geometry  of  the  cascade  had  been  determined  theoretically.  With  j 
the  aid  of  the  wake  traversing  method  the  total  head  losses  were  I 
measured,  the  discharge  angle  and  the  static  pressure  drop  across j 
the  cascade  were  checked,  and  this  both  with  and  without  a  chan-  ' 
nel  behind  the  cascade.  j 

i 

3  2 

j2)  Symbols  (see  figure  H2)  i 


3 


\ 

t 


« 


b 

c 

d 

r 

ji 

y 

po 

p 

q 

w 


CD 

cf 

po 1 ~po2 
*1 
Q 


Y 


Q 

P 


1-C 


blade  span  (=  300  mm)  { 

! 

blade  chord  (=  40  mm)  j 

distance  between  blades  (=  20  mm) 
radius  of  camber  line  of  the  circular  arc  sections 
axis  in  chordwise  direction 

axis  perpendicular  to  chordwise  direction  .  J 

total  head  j 

static  pressure  \ 

dynamic  head  J 

flow  velocity  j 

Reynolds  number  j 

t 

drag  coefficient  of  blade  jj 

a 

drag  coefficient  of  flat  plate  s: 

« 

i 

i 

loss  coefficient  of  cascade  I 


angle  between  flow  direction  ana  a  perpendicular  to  I 
the  cascade  axis  j 

stagger  angle,  angle  between  cascade  axis  and  a 
perpendicular  to  the  chord  j 

density  of  the  air  ] 

* 

dynamic  viscosity  of  the  air  ; 

efficiency,  of  the  cascade 
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Subscripts 


1 

upstream  of  blade  row 

2 

downstream  of  blade  row 

loc 

local  value 

S3)  Method  of 

test 

The  investigation  was  carried  out  with  the  small  cascade  tunnel 
of  the  Institute  of  Fluid  Mechanics  of  the  Technical  University 
Braunschweig  (director:  Professor  Dr.  H.  Schlichting ) .  A  view 
of  this  facility  is  given  in  Fig,H1 .  The  inlet  air  angle  could 
be  varied  by  turning  the  big  disk,  to  which  the  cascade  was  fixed. 

The  geometry  of  the  cascade  had  been  determined  theoretically 

with  the  aid  of  potential  flow  theory.  It  is  shown  in  .  The 

blades  had  circular  arc  shape  and  consisted  of  metal  with  a 

j  thickness  of  1,5  mm,  the  radius  of  the  camber  line  being  r=33,6mm. 
i'  | 

{  The  leading  edges  were  rounded  and  the  trailing  edges  sharpened  j 

thoroughly  to  get  good  airfoil  sections.  The  blades  were  set 

with  a  stagger  angle  of  y  =  2,5°.  The  total  deflection  of  the 

cascade  was  to  be  55°,  thus  the  inlet  air  angle  was  B.  =  25°,  ! 

I 

whereas  the  discharge  angle  was  expected  to  be  B^  =  30°.  Both 
the  inlet  and  discharge  directions  formed  an  angle  of  9°  with 
the  tangents  to  the  camber  line  at  the  leading  and  trailing 
edges  (see  Figure  H2) .  The  chord  length  was  c  =  40  mm,  the  dis¬ 
tance  between  the  blades  d  =  20  mm,  thus  the  solidity  d/c  =  0,5. 

>  ] 

The  cascade  consisted  of  20  vanes.  The  blade  span  was  b  =  300  mm. i 
Tost  was  performed  with  an  entering  velocity  of  w^  =  55  m/s 
corresponding  to  a  Reynolds  number  of  =  145  000. 

The  probe  used  was  a  conical  one  of  small  dimensions.  It  allowed 
simultaneous  measurement  of  total  head,  static  pressure  and  flow 
direction.  A  photograph  of  the  cascade  and  the  probe  is  shown  a 
in  Fig.  H3. 

First  of  all  the  cascade  was  investigated  alone.  It  was  set  with  j 
the  design  inlet  angle  of  B^  =  25°  ♦  The  total  head  loss  and  the  j 
discharge  angle  were  measured  behind  the  central  blade  at  a  I 

in  i 
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probe's  distance  of  one  chord  downstream  of  the  trailing  edge 
at  midspan  position.  They  were  taken  as  the  mean  values  over 
one  spacing  of  the  blades. 

j  After  this  a  channel  of  345  mm  height  and  a  length  of  its  center 
line  of  1,4  m  (see  FigureH2)  was  installed  downstream  of  the 
cascade  such  as  to  have  a  bend  with  a  total  deflection  of  55°. 
Measurements  were  taken  at  a  distance  from  the  trai ling .  edge  of 
x/l  =  6,0  in  the  middle  of  the  channel  ovc^  three  spacings.  They 
revealed  the  amazing  fact, that  the  losses  were  twice  as  high  as 
without  the  channel.  The  measurements  were  repeated  at  x/l  =  1,0 
to  have  a  better  possibility  of  comparing  with,  the  data  without 
channel.  Total  head  loss  was  only  reduced  slightly.  The  channel 
was  removed  again,  and  a  new  test  was  run  at  a  lower  Reynolds 
number  of  =  75  000,  since  it  was  presumed  that  the  Reynolds 
number  had  been  near  its  critical  value  with  the  tests  before, 
and  thus  a  slight  change  in  flow  pattern  caused  by  the  channel 
might  have  led  to  laminar  separation  of  the  upper  surfaces  (suc- 
liua  sides).  No  particular  rise  of  the  losses  was  found,  however. 
Bad  directional  setting  of  the  channel  could  not  be  the  reason 
for  the  high  losses  either,  since  it  had  been  done  carefully 
and  since  the  cascade's  discharge  angle  fitted  very  well  with 
the  channel's  direction. 

1 

j  An  answer  to  the  question  could  only  be  found  by  measuring  at 

| 

|  different  inlet  angles,  i.e.  by  determining  the  cascade's  polar 
l  curve.  Thus  additional  tests  were  run  with  the  isolated  cascade 
at  fi.  =  15°,  20°,  30°  and  35°. 

j  1 

r. 

>■ 

14)  Discussion  of  results 

5 

I. 

The  discharge  angle  fitted  very  well  with  theory,  it  was  found 
to  be  J39  =  29,5°,  i.e.  only  half  a  degree  less  than  was  to  be 

i 

j  expected  according  to  potential  flow  theory. 

!  The  total  head  losses  are  presented  in  Table  HI.  The  distribution 

l 

of  the  local  loss  coefficient 

\  r  =  ^Po  1~po2  Hoc 

j  Hoc  q1 
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'■'■i-:'*  channel,  as  measured  six  chords  behind  the  grid,.  is  shov.-r. 
in  Fir.H4  over  three  spaci,ngs.  The  wake  of  the  central '  blade  as 
measured  one  chord  downstream  of  the  cascade  is  shown  in  Fig,H5. 

As  can  be  seen,  the  losses  arc  C  =  0,18b  when  using  the  down- 
stream  channel,  but  only  C  =  0,092  without  it.  In  the  first  case, 
losses  reduce  slightly  to  C  -  0,166  when  measuring  only  one  chord 
behind  the  blade's  trailing  edge  instead  of  six  chords.  In  the 
second  case,  losses  increase  slightly  from  C  =  0,092  to  C  =  0,103 
when  lowering  Reynolds  number  from  R ^  =  1 4b  000  to  =  7b  000. 
Because  of  the  latter  result,  separation  at  the  upper  surface 
cannot  be  an  explanation  for  the  considerable  difference  between 
the  losses  of  the  cascade  without  and  with  the  channel. 

Variation  of  the  inlet  air  angle  showed  that  the  reason  for  the 
high  less  is  a  lower  surface  (pressure  side)  stall.  This  can  be 
seen  very  well  from  Fi&K5,  which  chows  the  wakes  at  different 
inlet  angles?1The  pressure  sides  of  the  wakes  get  larger  with 
decreasing  inlet  angle,  the  tosses  thus  increase. 

A  better  understanding  of  the  unaer tying  flow  mechanism  yields 
the  polar  curve  of  FigHb.  The  total  head  loss  with  channel  cor¬ 
responds  to  the  loss  of  a  cascade  without  channel  that  has  a 
7,3  smaller  inlet  angle,  namely  =  17,7  .  Obviously  installing 
a  channel  behind  the  cascade  causes  rather  a  severe  change  of 
its  inlet  angle .  Thus  the  idea  that  a  channel  downstream  of  a 
cascade  does  not  change  the  flow  conditions  in  front  of  the 
cascade,  if  only  the  channel  slopes  back  with  its  discharge  angle j 
does  not  appear  to  be  correct.  Maybe  it  is  with  blades  of  much  ! 
larger  chord.  But  in  the  present  case  one  must  adopt  another 
concept,  namely,  first  of  all,  imagine  a  bent  channel  without 
any  grid.  The  streamlines  then  will  be  curved  in  some  particular 
way.  When  inserting  a  cascade,  this  will  experience  a  different 
inlet  air  angle  than  if  it  were  isolated,  according  to  the  al¬ 
ready  existing  curvature  of  the  streamlines.  Thus  the  change  of 
the  grid's  inlet  angle,  when  using  a  downstream  channel,  seems 
to  be  quite  reasonable.  This  fact  may  or  may  not  change  the  loss 
across  the  cascade,  depending  on  its  geometry.  In  the  present 
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case  it  led  to  a  big  change,  since  the  inlet  conditions  of  the 
grid  were  bad  already  at-  the  design  angle  of  fi|  =  21°.  Namely, 
the  polar  curve  (Pig.H6)  shows  that  the  blades’  lower  surfaces 
were  stalled  already  then.  The  reason  for  this  is  the  fact  that 
the  flow  direction  forms  an  angle  of  9°  witn  the  tar  ■.•at'  to  the 
camber  line  at  the  leading  edge  (see  ?ig.H2),  thus  th  stagnation 
point  being  located  at  the  blade's  upper  surface.  Even  a  small 
decrease  of  inlet  angle  will  lead  to  a  considerable  increase  in 
total  head  loss  then. 

2ig.H6  shows  also  that  the  cascade  works  rather  far  away  from 
its  optimum  configuration  which  exists  with  an  inlet  angle  of 
8-j  =  32°  instead  of  13^  =  23°.  The  loss  coefficient  then  is  only  ! 
C  =  0,045  instead  of  C  =  0,092.  It  does  not  appear  to  be  very 
favourable,  ho  never,  to  have  the  isolated  grid  work.  at  its' 
optimum  inlet  angle.  Namely,  as  was  shown  by  the  experiments, 
the  set-up  with  the  channel  then  would  work  a z  a  smaller  inlet 
angle,  i.e.  in  a  region  of  higher  losses,  minimum  loss  with  the  j 
channel  configuration  will  probably  be  attained  at  inlet  angles  | 
n the  isolateu  grid  that  are  5°  to  7°  higher  than  its  optimum  ! 
angle,  i.e.  at  about  3^  =  38°. 

Fig ,H7  shows  the  efficiency  of  the  cascade  plotted  versus  the  ! 

i 

inlet  angle.  It  is  defined  by  j 


*2  '*  1  "  CD  *  d 

In  the  case  of  a  cascade  with  zero  stagger  (see  Ref.  HI) 


where 


i  -  c 


po1'Po2 


The  graph  of  Fig.H7was  calculated  with  equ.(3)  since  the  stagger 
of  only  y  =  2,5°  can  be  neglected. 

?ig.H7  shows  also  the  efficiency  of  a  cascade  of  flat  plates  with 
zero  deflection  anu  when  separation  is  not  taken  into  account. 

It  was  found  by  introducing 


a 


Insiitut  fur 
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8km  Nr.  7 


CD  =  cf  ~  0  >  00 6 6  , 

the  value  for  turbulent  flow  at  a  Reynolds  number  of  =  190  000. 
The  flat  plate  efficiency  is 

1,  flat  =  0,987  '  » 
plate 

whereas  the  cascade's  optimum  efficiency  is  y  =  0 ,  Qpd ,  and  taat 
one  at  the  design  inlet  angle  is  only 

y  =  0,903 


without  a  channel  behind  the  cascade. 

Suggestions  for  improvement 

It  is  of  no  use,  however,  to  increase  the  cascade's  inlet  angle 
since  the  total  deflection  wpuld  be  greater  than  99°  in  that  case. 
Improvement  can  only  be  obtained  by  changing  the  cascade's  geo¬ 
metry.  The  trailing  edge  geometry  should  be  kept  the  same  because 
the  discharge  angle  fitted  well  with  theory.  But  the  leading 
edge  geometry  should  be  altered  such  that  the  flow  direction  and  j 
the  tangent  to  the  camber  line  do  not  form  an  angle  of  +9°  as 
They  do  with  the  present  configuration,  but  an  angle  of  -5°  to  '  \ 
-7°.  The  set-up  with  the  channel  will  work  then  just  at  the  - 

cascade's  optimum.  j 


t. 9 )  Conclusions 

i  The  loss  coefficient  of  a  cascade  of  circular  arc  airfoils  was 
! 

investigated  experimentally  with  and  without  a  channel  downstream 
I  *  j 

|  of  the  cascade.  The  turning  angle  was  checked.  It  was  found  that 
the  discharge  angle  compared  very  well  with  potential  flow 
j  theory.  But' the  losses  were  rather  high  since  the  cascade  did 

t 

j  not  work  at  its  optimum  configuration.  The  set-up  with  the 

j  channel  yielded  even  higher  losses  since  the  inlet  air  angle  of 

the  cascade  was  shifted  by  the  channel  in  an  unfavourable  way. 

I 

|  The  geometry  of  the  cascade  should  be  altered  by  leaving  the 
trailing  edge  configuration,  but  diminishing  the  camber  such  as 
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to  have  an  ungle  of  -5°  to  -7°  between  the  flow  direction  and 
the  tangent  to '  the  camber  line  at  the  Leading  edge,  instead  of 
+S°  as  with  the  present  vanes. 
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Resultant  Loss  coefficients 


\ 

I  B.'olf  Nr. 


Table  h  1 


Pi 

[cleg] 

without  channet 
x/c  =1,0 

l"'"  '  ■  _ 

with  channet 

Rj  =150000 

*7 

£ 

x/c 

15 

195000 

0,192 

— 

20 

195000 

0, 102 

- - 

- - 

25 

195000 

0,092 

1,0 

0,166 

75000 

0, 105 

6,0 

0,195  V 

JO 

195000 

0,052 

- - 

55 

195 000 

0,053 

- — —  Ik 

— 

- - 

nlThis  imlue  is  the  mean  Loss  coefficient  over  three  spacings 
y/hiLe  all  the  others  were  obtained  o/er  one  spacing  only  ' 


U 


x/c:  Distance  of  probe  behind  trailing  edge  in 
chordwise  direction 


125 


SECTION  I 


COMMENTS  ON  CASCADE  TEST  RESULTS 


The  tests  made  by  the  "Institut  fur  Stromungsmechanik  M  at 
Braunschweig,  Germany  are  of  great  value  to  determine  the  optimum 
cascade  for  use  in  VTO.  It  is  to  be  stated  categorically  that  a 
propeller  or  fan,  under  no  circumstances,  should  be  exposed  to  a 
non- axial  stream.  Not  only  are  the  mechanical  stresses  excessive 
but  also  the  efficiency  of  the  propeller  is  drastically  reduced  under 
the  very  condition  for  which  the  highest  efficiency  is  required. 
Propellers  in  ducts  have  been  tested  adequately  for  any  and  all 
purposes.  Nevertheless,  there  is,  of  course,  still  need  for  ordinary 
model  testing  as  basis  for  actual  design  and  construction  of  VTOL 
airplanes,  in  particular  to  insure  compatability  of  the  propeller  and 
duct. 


With  propeller-duct  and  ordinary  airplane  testing  quite  fully 
explored  it  was  considered  desirable  to  obtain  design  information  on 
cascades  as  the  only  field  not  adequately  explored.  It  was  for  this 
reason  that  the  present  test  series  was  conducted. 


I 


The  experimental  cascade  had  been  designed  from  exact  theory 
(Section  G  of  this  volume).  The  test  data  revealed  the  interesting  fact 
that  the  duct  walls  should  have  been  eliminated.  In  other  words,  the 
cascade  should  have  been  tested  as  a  free  stream  deflector  or,  if  in  a 
duct,  the  duct  walls  should  have  been  streamlined.  Fortunately,  because 
of  the  extensive  data,  including  wake  surveys,  it  is  quite  simple  to 
obtain  the  desired  information  on  the  efficiency  of  properly  designed 
baffles.  The  actual  result  may  be  obtained  as  follows: 


The  channel  on  the  inlet  side  which  was  present  in  all  the  tests 
forced  a  straightening  of  the  flow  (parallel  to  the  walls)  so  that  the  flow 
enters  the  cascade  near  zero  degrees  instead  of  the  theoretically  pre¬ 
scribed  9°  (see  Figure  H2  in  Institute  report).  As  the  flow  is  prevented 
from  entering  at  the  proper  angle  it  stalls  on  the  pressure  side.  In  fact, 
the  highest  efficiency  is  reached  when  the  inlet  angle  ^  =  approximately 
32°  (see  Figure  H7)  which  gives  an  angle  of  attack  of  the  air  with  the 
leading  edg^  of  the  blades  of  only  2°.  Notice  the  drastic  reduction  in  the 
loss  on  the  pressure  side  (Figure  H5  )  at  /&  =  30  .  The  efficiency  reached 


about  95%  in  spite  of  the  fact  that  one  has  now  a  diverging  channel.  The 

cross  section  of  the  duct  on  the  outlet  side  is  proportional  to  cos  25°  and 

o  •  866 

on  the  inlet  side  to  cos  30  or  the  ratio  is  -  which  corresponds  to 

.94 

more  than  7%  divergence  instead  of  the  originally  intended  convergence  of  7% 


.  o 

with  the  outlet  angle  at  A  a  30  and  the  inlet  at  jS  -  25  (in  reverse  order), 

'2  1 
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It  is  proper  to  expect  the  test  conforming  with  the  theoretically 

prescribed  inlet  and  outlet  angles  to  give  proper  condition  for  the 

suction  side  (FigureH5 curve  for  \  *  25 °  without  outlet  channel). 

Further,  it  is  proper  to  use  the  curve  for  ^  ^  -  30°  for  the  pressure 

side  since  this  side  is  not  too  much  affected  by  the  divergence  of  the 

channel.  We  obtain  thus,  from  FigureH5by  drawing  an  internal  contour 

curve,  a  rather  accurate  measure  as  to  the  less  in  the  ideal  case.  This 

loss  seems  to  be  between  1/2  and  2/3  of  that  of  the  best  case  observed. 

o 

As  the  lowest  loss  measured  (in  the  diverging  channel  with  a  30  ) 
was  4.  5%  the  ideal  loss  is  thus  estimated  as 

-  2.25  -  3.% 

This  compares  with  a  skin  friction  loss  of  1.  3%  as  the  absolute  minimum. 

We  may,  therefore,  state  in  final  conclusion  that  carefully 
designed  baffles  may  reach  an  efficiency  of  97  to  98%.  Thus  such 
elements  may  be  used  freely  for  any  purpose  of  directing  the  flow  in  a 
VTO  airplane.  The  incurred  loss  is  much  less  than  that  caused  by  a  non- 
axial  inflow  into  a  propeller. 

It  is  of  interest  to  note  that  Dr.  Schlichting  had  given  the  value 
of  98%  as  the  expected  value  for  the  efficiency  in  the  initial  discussion 
of  the  contract  negotiations  on  the  cascade  tests  at  Gottingen. 
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